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Introduction 



The birth and the development of quantum mechanics is a remarkable example 
of the scientific method introduced in the 17th century by Galileo Galilei. He 
defined this method as the composition of two fundamental steps: "sensate 
esperienze" (sensory experience) and "necessarie dimostrazioni" (necessary 
demonstrations) PQ . 

Unexplained phenomena like photoelectric emission, black body radiation, 
atomic emission lines, etc., encouraged the definition of a new theory called 
quantum mechanics. This is the first step. However, simultaneously with the 
birth of the theory, scientists started to ask what are all the possible impli- 
cations of the new quantum mechanical laws. This is the second step. They 
discovered that these implications are very interesting and unusual, e.g. a 
system can be in a superposition of states, the result of a measurement can be 
predicted only as a probability, position and momentum of a particle cannot be 
simultaneously measured with infinite precision, the effect of a measurement 
on a system can instantaneously affect the state of another spacelike separat- 
ed system, etc.. Moreover, during the last century, other interesting questions 
were posed: "Is quantum mechanics a complete theory? Or is it a manifes- 
tation of an underlying hidden variable structure?", "Can we use quantum 
mechanics to transmit information or to perform computations, with better 
performances with respect to classical methods?", "If quantum mechanical 
laws are true for microscopic systems, why should not be true also in the 
macroscopic world?", etc.. This thesis is focused on the second step of the 
scientific method which Galilei called "necessarie dimostrazioni". In fact we 
consider some important implications and predictions of quantum mechanics 
and we study how they can be experimentally verified. 

In particular, in the first part of the thesis, we consider the quantum 
mechanical phenomenon of entanglement, according to which the state of a 
system is so much correlated to the state of another spacelike separated system, 
that an operation performed on the first one causes a distant action on the 
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second one. 

Within the formalism of quantum optics, we introduce also the technique of 
quantum teleportation. This is a smart way of using entanglement and classical 
communication, to transmit an exact copy of a quantum state between two 
spacelike separated interlocutors (usually called Alice and Bob) . We determine 
also which are the best local operations, that Alice and Bob can perform at 
their sites in order to optimize the teleportation of a coherent state. For a 
fixed amount of entanglement (negativity), we find upper and lower bounds for 
the maximum of the teleportation fidelity and we prove that this maximum is 
reached iff the correlation matrix of the shared entangled state is in a particular 
normal form. 

Usually it is said that quantum mechanics is the study of physical systems 
whose dimensions are very small such as atoms, electrons, subatomic particles, 
photons, etc.. However this is not completely true, since quantum mechanics 
is a new way of looking nature and it is a theory which describes an electron 
as well as a planet. One of the aim of this work is indeed to give theoreti- 
cal predictions of quantum effects in macroscopic systems. In particular, in 
the second part of the thesis, we will investigate optically driven mechanical 
oscillators: a new emerging research field named quantum optomechanics. 

We consider an optical cavity in which one of the two mirrors is free to 
oscillate around its equilibrium position. The motion of the mirror is coupled 
with the cavity light field by radiation pressure. We will show that this cou- 
pling can optically cool a vibrational mode of the mirror and can entangle the 
state of the macroscopic mirror with the intracavity field. Moreover, we show 
that appropriately filtering the cavity output field, different optical modes can 
be extracted. These modes possess a significant amount of entanglement with 
the mirror, which can be also higher than the entanglement of the intracavity 
field. 

In the last chapter we consider also a cavity driven by two different lasers. 
We show that, by choosing opposite detunings and equal laser powers, a stable 
regime can be reached. In this regime, optomechanical entanglement between 
the laser fields and the mirror is improved. Also the two optical output fields, 
if appropriately filtered, are robustly entangled even at room temperature. 

Finally we apply the classification criterion of [53] to show that, both in 
the monochromatic and in the bichromatic setup, we can find two optical 
output sidebands which together with the mirror vibrational mode are in a 
fully tripartite-entangled state. 
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CONTINUOUS VARIABLE 
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INFORMATION 



Capitolo 1 

Quantum mechanics 



It is well known from quantum mechanics, that a pure state of a quantum 
system can be represented by a ray of a Hilbert space associated to the system. 
This ray is written in the Dirac notation [6] with a ket vector \ip) and it is said 
to be pure since the uncertainty that it possesses is due only to the Heisenberg 
principle and not to our lack of information about it. However, in general and 
especially in quantum optics, we do not know exactly which is the state vector 
of a system, but we know only the probability P{a) that the system is in the 
state \i(i a ). In this case the system is said to be in a mixed state. For this 
reason it is better to use a mathematical object, firstly introduced by Von 
Neumann [S], which is more powerful than the usual ket vector: the density 
operator. 

1.1 Density operator 

Definition 1. We associate to every state a density operator p defined through 
the following convex sum of projectors 



where P(a) is the probability that the system is in the pure state vector \i/j a ) 
and therefore we must require ^2 a P(a) = 1. 

This new formalism associates to each possible state of a system an op- 
erator p which contains all the quantum and statistical informations that we 
know about the system. 

The expected value of a generic observable O is given by 




(1.1) 



a 




(1.2) 



a 
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and we can write it in a compact form as 

(0) = Tr{pd]Q. (1.3) 
Three important properties follow form the definition, 

P = p\ (1-4) 
Tr{p} = 1, (1.5) 
P>0. (1.6) 

Definition 2. For every state p we can define [5] the Von Neumann entropy 
H(p) as 

H(p) = -Tr{plnp}. (1.7) 

The Von Neumann entropy can be considered as a measure of the disorder 
of a given quantum state. Moreover, at thermal equilibrium, it coincides with 
the classical thermodynamic entropy S of the system (up to a multiplicative 
factor given by the Boltzmann constant Kb)- 

1.2 Pure and mixed States 

Definition 3. Given a density operator p = Y2 a -^WIVO {ipa\, if the set of 
vectors {(V'a)} is made of only one element, the state is pure, otherwise it is 
mixed. 

Theorem 1. A state p is pure iff p 2 = p. 

Dimostrazione. p 2 = = 9- 

<^= yo is Hermitian and therefore it can be diagonalized 



J2*n\n)(n\. (1.8) 



The hypothesis p 2 = p implies A n = 0, ±1, but for the normalization condition 
Tr{p} = 1 and p > the only possible eigenvalues are A n = 5 n ^ and therefore 
p=\k)(k\. □ 

Another way to distinguish pure states from mixed states is provided by 
the trace of p 2 . 

Theorem 2. A state is pure iffTr{p 2 } = 1, while is mixed iffTr{p 2 } < 1. 



1 The trace of an operator A is defined as Tr{A} = (n\A\n), where {|n)} is an arbitrary 
orthonormal basis of the Hilbert space. 
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Dimostrazione. 

Tr{p 2 } = Ys P » P b\^aUb)\ 2 <Y, P - P b = !. (1-9) 

a,b a,b 

and the inequality becomes an equality only when P a Pb = #a,&, that is if 

p = |V>><# □ 

Finally one could prove also the following theorem, showing that the purity 
of a state is connected with the amount of information that we have about it. 

Theorem 3. A state p is pure iff H(p) = 0, where H(p) is the Von Neumann 
entropy defined in fli. 7| ). 

1.3 Dynamics 

In the Schrddinger picture the time evolution of a state is determined by the 
Hamiltonian H of the system through the Von Neumann equation (the density 
matrix equivalent of the Schrodinger equation), 

>) = al*.fl d-io) 

which can be formally integrated giving 

p{t) = exp(-^)p(0) exp^Ht), (1.11) 

while all the observables do not depend on time. 

A completely equivalent formulation of the system dynamics is given by 
the Heisenberg picture, in which the state p remains constant in time while 
the observables evolve obeying the Heisenberg equation 

which integrated gives 

A(t) = exp(lHt)A(0) eM-^Ht). (1.13) 
n n 

Finally we can give also another formulation which is intermediate between 
the two that we have seen: the interaction picture. Let us suppose that the 
Hamiltonian (in the Schrodinger picture) can be written as H = Hq + Hi, 
where usually Hi is an interaction/perturbation term. We want to stay in 
a reference frame such that the effect of Hq cancelled. This can be done by 
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making the following transformation to the operators and to the interaction 
Hamiltonian H±, 

Aj(t) = exp(-^H t)A(t)exp(^H t), (1.14) 
tfj(t) = exp(-^ot)£iexp(^ i). (1-15) 
If we derive (|1.14|) . using (|1,12|) . we get 

±A I (t) = ^[A I (t),H I (t)], (1.16) 

which is very similar to (|1.12|) . but in this case the evolution depends only on 
the interaction Hamiltonian Hj. 



Capitolo 2 

CV systems in phase space 



A quantum system described by observables with a continuous spectrum is 
said to be a continuous variable (CV) system. For example, the state of a 
particle in a generic potential is a CV system, since it is described in terms of 
its position and momentum which are continuous observables. 

All the quantum mechanical features of a system are very fragile if this 
system interacts with the environment, since the effect of decoherence destroys 
these quantum effects in a very short time. For this reason, thanks to its 
limited interaction with the environment, the electromagnetic radiation is a 
good candidate for observing quantum phenomena. It is well known that, in 
the canonical quantization of the electromagnetic field, a radiation mode is 
formally equivalent to a harmonic oscillator. Therefore when we talk about a 
CV system we usually refer to one or more harmonic oscillators described (if 
they are isolated) by the simple Hamiltonian 

1 N 

H = -J2^k(f k + ql), (2.1) 
k 

where q^ and pk are dimensionless position and momentum operators satisfy- 
ing the canonical commutation relations (CCR): [qkiPi] = i&u and [qk,qi] = 
\pk,Pi] = 0. If we define the quadratures vector R = (qi,pi, q2,p2, ■ ■■(1n,Pn) T 
we can express the CCR in matrix form 

[A fc ,i?j]=in w , (2.2) 

where Q is the 2N x 2N symplectic matrix given by 
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Since we deal with harmonic oscillators, it is useful to introduce also the 
annihilation and creation operators which are respectively 

Ofc = -j= , (2.4) 



a fe = (2.5) 



and satisfy the bosonic commutation relations [dfc,d[] = 5u and [a^, d;] = 

[4,af] = 0. 

If we consider a single oscillator, it can be shown that the ground state 
of the Hamiltonian satisfies d|0) = and that all the other eigenstates (Fock 
states) are given by 

|n) = -=|0>. (2.6) 
Vu! 



Fock states have relatively large variances on the position and momentum op- 
erators. A quantum description which is close to the classical one is given by 
the coherent states which are characterized by equal and minimum uncertain- 
ties in the quadratures Aq = Ap = These states are the eigenstates of the 
annihilation operator a\a) = a\a) or equivalently \a) = exp(adt — a*d)|0), 
where a is the mean amplitude (d) of the state. 

The states with minimum uncertainties but different variances in q and p 
are the squeezed states \a,r) = exp(ad^ — a* a) exp(r 2 d^ 2 — r*d 2 )|0). These 
states (coherent and squeezed) are very important in quantum optics, however 
we do not give a deep treatment here since they are particular cases of the 
most general class of Gaussian states studied in the next section. 

If we want to give a quantum description of a system of oscillators we need 
to work with operators acting on infinite Hilbert spaces and this can create 
some problems. To overcome these problems it is better to use an equivalent 
representation in which we use c-number functions instead of operators and 
we work with a finite number of dimensions: the phase space. 

In classical mechanics, the state of a particle, or an ensemble of particles, 
is represented by a probability density P(q,p) in phase space. With this 
distribution is possible to calculate the mean value of every observable, 

(0(q,p)) = J P(q,p)0(q,p)dxdp, (2.7) 

in fact, if we power expand the observable 0(q,p), eq. f|2.Tj) reduces to a sum 
of mean values of terms like q n p m 

(0(q,p)) = J> nm I' q n p m P(q,p)dqdp. (2.8) 

n,m 



Chapter 2 - CV systems in phase space 



9 



If we come back to the quantum world, it is natural to try to define a dis- 
tribution W(q,p) in phase space similar to the classical one P(q,p). However 
if we deal with quantum states two problems arise. The first is that we have 
to consider the intrinsic quantum uncertainty of a state, so that for example a 
distribution which is too much localized is forbidden by the Heisenberg princi- 
ple. Another problem is that in (|2.8p the variables q and p commute while the 
quantum operators q and p no not commute. For this reason we have different 
quantum distributions depending on the order of the variables products that 
we want to calculate. 



2.1 The Wigner function 

The fundamental operator that we need in the passage from the Hilbert space 
to the phase space is the Weyl operator defined as 

Wt = e ienA , (2.9) 



where f2 is the symplectic matrix, R is the quadrature vector and £ £ R . It 
can be shown that the CCR ([2.2D are equivalent to the Weyl relations 



W^W 62 = e^'Wfr+fc. (2.10) 

As we can describe a CV system using the quadrature operators satisfying the 
CCR, in the same way we can give an equivalent description in terms of Weyl 
operators satisfying (|2.1U|) . Using these operators we can define something 
which is similar to a Fourier transform but it is applied to operators instead 
of functions: the Fourier-Weyl Transform. 

Definition 4. The Fourier-Weyl Transform is a map from 7i to £ 2 (R 2Ar ) 
which associate to each operator A the function 

A(£) = Tr(AWz). (2.11) 

The inversion of (|2.11|) give us the possibility to express any operator as a 
linear combination of Weyl operators, in fact we have 

A =J^wJ MOW^dC (2.12) 

In analogy with the classical one, it can be shown also the quantum Parseval 
theorem 

Tr(iB) = J A(0*B(0(%. (2.13) 
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The most important application of the Fourier- Weyl Transform is when 
we use it to transform the density operator. In this case the result is called 

characteristic function 

X (0 = Tr(p^), (2.14) 

because of the analogy with the classical characteristic function which is the 
Fourier transform of a density distribution. For the Parseval theorem (|2.13p . 
the properties of the density operator that we have seen in the previous chapter 
imply 

X(0) = 1, (2.15) 
|X(£)| 2 # < 1, (2-16) 



N 



(27T) 

where the last integral is equal to 1 iff the state p is pure. 

Finally, an important property of the characteristic function is that, through 
its derivatives, we can express all the moments of the density operator. In fact, 
if we define R = QR, we have 



( _j\".+m ara+m 



(2.17) 



where {,} sy m means the symmetrized product. 

So, thanks to the characteristic function we can map a density operator to 
an element of C 2 , but if we want something which is analogous to a probability 
distribution in a phase space we need the inverse Fourier transform of x(v) 

W(0 = J xOrte-^de, (2.18) 

which is called the Wigner function. This function is very similar to a proba- 
bility distribution but it has some different properties, for example it can take 
also negative values. However it is true that, like a density probability, the 
Wigner function can be used to calculate the mean values of symmetrically 
ordered operators 

Tv[p{(R' k r, (RTUm] = J mrmm- (2.19) 

2.2 Gaussian states and Gaussian operations 



Definition 5. A quantum state p associated to N harmonic oscillators is 
Gaussian if its characteristic function is a Gaussian. 
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If we define the displacement vector d = Tr{pR}, and the symmetrically or- 
dered covariance matrix (CM) V/y = Tr{/3[i?£.— df., Ri — di]+}, the characteristic 
function of a Gaussian state is 

where we can see that the first and the second moments completely define the 
state. Prom the definition (|2.18p . the associated Wigner function is 

= ^| F |-i e -(€-«0 T n T v- l n«-«O j ( 2 . 2 1) 

where |V| is the determinant of the CM. 

We must observe that not every Gaussian function can be a Wigner func- 
tion of a physical state. In fact the Heisenberg principle imposes a constraint 
on the correlation matrix [9]. 

Theorem 4. A correlation matrix V correspond to a physical Gaussian state 
iff 

V + in>0. (2.22) 

We have seen that through the Wigner function we can have a rigorous 
phase space representation of the state p of a system. Now we may ask: if 
we make a physical operation on the state p what is the effect on the Wigner 
function? We will answer to this question considering only the set of Gaussian 
operations, which are the transformations p — ► e{p) which maps Gaussian 
states into Gaussian states. Since a Gaussian state is determined by its first 
and second moments we can completely define a Gaussian operation describing 
its effect on the correlation matrix V and on the displacement vector d. 

2.2.1 Displacement transformations 

Definition 6. A displacement transformation is a translation in phase space 
parametrized by the displacement vector r S M. 2N . It changes only the first 
moments of a Gaussian state 

d — ► d' = d + r. (2.23) 

The corresponding operator in the Hilbert space is the Weyl operator, we 
have indeed that 

W T RWl =R + t. (2.24) 

Any Gaussian operation can be seen as a composition of a transformation 
leaving the first moments unchanged followed by a displacement. For this 
reason, all the next operations we are going to study will act only on the 
correlation matrix V of the Gaussian state. 
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2.2.2 Symplectic transformations 

Definition 7. A symplectic transformation is parametrized by a 2N x 2N 
matrix S which transforms the set of coordinates R to R! = SR, leaving the 
CCR unchanged 

sns T = n. (2.25) 

The effect on the CM of a Gaussian state is 

V — >V = SVS T . (2.26) 

The set of these transformations forms the real symplectic group Sp(2N, R), 
which is a subgroup of the real special unitary group SU(2N, R) since we always 
have det S = 1. 

Theorem 5 (Metaplectic Group). Given a symplectic transformation S, we 
can always find a unitary operation Us in the underlying Hilbert space which 
realizes S in the phase space. In other words, for every S £ Sp(2N, R) there 
exists an (up to a phase unique) unitary Us such that 

UsWfil = W S £, V£ e R 2N . (2.27) 

The set of all these unitaries Us forms the metaplectic group Mp(2N), 
which is a two-fold covering of Sp(2N, R). In fact, given a symplectic matrix 
S we always have a ± ambiguity in the choice of Us- 

Theorem 6 (Euler Decomposition). Every symplectic transformation S £ 
Sp(2N, R) can be decomposed in the following way 

S = R 2 U$dAri, (2.28) 

where Dk are 2x2 symplectic diagonal matrices Dj~ = diaglr^,^ 1 }, while 
R12 S Sp(2N, R) n SO(2N, IR) are symplectic rotation matrices. 

The unitary transformations Ur 12 £ Mp(2N) associated to R12, can be 
physically realized with passive optical elements (beam splitters and phase 
plates). The matrices instead correspond to single mode squeezing opera- 
tions, which must be realized using active optical elements (non linear crystals, 
homodine measurements, etc.). 

Theorem 7 (Williamson [TU] ) . Every positive definite 2N x 2 N matrix V can 
be diagonalized by a symplectic transformation S € Sp(2N, R) in the following 
form 

SVS T = diag[si, si, s 2 , s 2 , ....sn,sn]- (2.29) 
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The elements s& are called the symplectic eigenvalues and from 112. 22\) we can 
say that V is a CM of a Gaussian state iff s^ > 1 for every k. 

2.2.3 Trace preserving Gaussian CP maps: Gaussian channels 

Definition 8. A trace preserving Gaussian completely positive map (TGCP) 
is parametrized by a pair of 2N x 2N matrices S and G = G T > 0, which 
satisfy the constraint G+i£l — iSQS T > 0. The effect on the CM of a Gaussian 



Physical Gaussian operations like amplification/attenuation or any other 
one in which we make no measurements is a TGCP map. These map are also 
called Gaussian channels since for example a dissipative medium in which 
there is an interaction with other external modes is modelled by a TGCP 
map. 

Finally, another definition of a TGCP map comes from the dilation theo- 
rem, which asserts that every Gaussian channel can be seen as the effect on a 
subsystem (partial trace) of a global symplectic map acting on a larger Hilbert 
space. 

2.2.4 Non-trace preserving Gaussian CP maps 

The most general physical operation acts on the density operator as a com- 
pletely positive map p — > e(p). To characterize a CP map it is useful to apply 
the Jamiolkowski isomorphism which establishes a one to one correspondence 
between CP maps e E H and positive definite density operators in the doubled 
Hilbert space p £ ^H®H [53] . 

Given the CP map we can construct a physical state p £ starting from a 
two-mode infinitely squeezed state 



state is 



V 



V = SVS T + G. 




oo 



\EPR) = lim cosh(r) \ tanh(r) n |n) (g> \n) 



(2.31) 



r— >oo 



n=0 



and then applying e to the first mode 



p £ = (e®I)\EPR). 



(2.32) 



The inversion of such correspondence is also possible. In fact, given a state p £ 
it can be shown that, if we use it as a channel for the teleportation of a state 
p, we can get (probabilistically) at the output e(p). If the CP map is Gaussian 



14 



Chapter 2 - CV systems in phase space 



then also the state p e is a Gaussian state and its CM completely characterizes 
the map. 

Definition 9. A Gaussian completely positive map (GCP) is parametrized by 
a 4N x 4N correlation matrix V corresponding to a physical state 



r 



ri r 12 

Tl2 ^2 



> -in e o. 



The effect of the map on the CM of a Gaussian state is 

1 



V 



V 



r i2 A 



Ar 2 A + v 



AIT 2j 



(2.33) 



(2.34) 



where A = diag[l, — 1, 1, —1, . . . ]. 



Any physical Gaussian operation, including also measurements and clas- 
sical communication, is a GCP map and can be described by (|2.34p . For 
example given a bipartite state with CM 



V 



A C 
C T B 



(2.35) 



then the projection of the second mode B into the vacuum |0)£0| is associated 
by the Jamiolkowski isomorphism to the physical state of C 

Ti = 



1,2 



r. 



cosh r \ 

coshr J 

sinhr 

— sinh r 

/ coshr \ 

coshr 

10 

1 



(2.36) 
(2.37) 

(2.38) 



where r — ► oo. From (|2.34j) . it can be shown [HI] that the action of the map 
on V is 



V 



V = A-C- 



1 



B + I 



:C T . 



(2.39) 



Finally, we remark that the formula given in (|2.30p . defining a TGCP map, 
can be obtained particular limit of (|2.34p [34J. 



1 We assume to loose the measured mode, therefore the dimension of T is not (4 + 4) x 
(4 + 4), but (2 + 4) x (2 + 4). 



Capitolo 3 

Entanglement 



Erwin Schrodinger, in his famous publication about the Schrodinger's cat para- 
dox of 1935, firstly introduced the word entanglement to define a quantum 
correlation that can exists between two or more systems, so that the state of a 
system is affected by the state of the others also if these systems are spatially 
separated. Entanglement is a phenomenon which is expected form the pos- 
tulates of quantum mechanics and it is based on the quantum superposition 
principle. In classical mechanics, superposition is not expected and therefore 
entanglement cannot exist, for this reason it seems so strange and paradoxical 
to our daily vision of the world. 

3.1 Historical and mathematical points of view 

The problem of entanglement, also if this word did not exist yet, was pointed 
out by Einstein, Podolsky and Rosen through the EPR paradox [3], which 
seemed to bring into question the foundations of quantum mechanics. 

Let us take a look to it. Einstein et al. considered two particles such that 
x\ — X2 = and p\ +P2 = 0. This is not forbidden by the uncertainty principle, 
since x\ — xi = X- commutes with f>\ + p2 = P+- So, it is possible to realize 
a physical state, the EPR pair, such that the previous operators are equal to 
zero with infinite precision: AZ_ = AP + = 0. 

If we postulate the hypothesis of locality, according to which a measure 
made on a particle cannot affect the state of another space-like separated 
particle, then the paradox arises. In fact, by measuring x\ and p2, it seems 
possible to know with arbitrary precision also x<i and p\, but this would violate 
the Heisenberg principle. The paradox can be avoided if we allow quantum 
mechanics to be a nonlocal theory (as it has been actually demonstrated), so 
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that a measurement on a particle changes the state of the other one and the 
previous argument cannot be used. 

Even if the aim was to criticize the completeness of the quantum theory, 
Einstein at al. actually centered the heart of quantum mechanics: entangle- 
ment. The principal feature of entanglement is indeed the possibility to have 
nonlocal effects and the EPR pair is still now the cornerstone of every CV 
information theory. 

We have briefly seen the historical origin of entanglement, now we give a 
rigorous mathematical definition of it. 

Definition 10. Given N subsystems 1,2, ... N, with the associated Hilbert 
spaces H\,Hi, ■ ■ ■ Hn, according to the postulates of quantum mechanics, the 
whole system is an element of the tensor product of the iV Hilbert spaces. A 
state p £ TL\ ® Ti.2, ■ ■ ■ ® Hn is separable if it can be expressed as a convex 
sum of operators like (g) , ■ ■ ■ ® , where 6 Hi . That is 

P = J2PkPk ] ®Pk ) >---®Pk N) > Pk>^ J>fc = l- (3.1) 

k k 

Otherwise the state is said to be entangled. 

We observe that, given a composite state, different but equivalent decom- 
positions in terms of £ Hi exist. For this reason, if a state is not written 
explicitly in the form (|3.1|) . we cannot be sure that it is entangled, in fact 
there could be a different decomposition showing that the state is separable. 

For Gaussian states we can give an analogous definition which depends on 
correlation matrices instead of density operators. 

Definition 11. A Gaussian state p G Hi <8> H2, • • • Hjv whose CM is V is 
separable if a set of Gaussian states g Hi with respective CM denoted by 
V 1 exists, such that 

V >V {1) @V {2) ■■■®V^ N) . (3.2) 
Otherwise the state is entangled. 

3.2 Separability criteria 

The previous formal definitions (|3.1)3.2p are very simple but they are often 
useless from a practical point of view. In fact, we usually need an opera- 
tional criterion to distinguish separable from entangled states. In the following 
subsections we give some of the most relevant separability criteria known in 
literature. 
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In the next chapters we will consider almost always lxl modes bipar- 
tite state, therefore the descriptions of the following criteria will be given for 
bipartite states. However, these descriptions could be easily extended to N 
parties states without introducing any conceptual difference. 

3.2.1 Peres-Horodecki criterion 

If p is a density operator of a physical state, then given an arbitrary orthonor- 
mal basis we can write the state as p = ^2ij fq\i){j\, where the Hermitian 
matrix fp- is the representation of p in the given basis. 

We define the transposition operation p — ► p T as the map which associate 
each density operator represented by p? to the operator represented by p*-. 

We observe that if p is a physical state than also p T is a density operator 
of a physical state, in fact it satisfies the same three conditions given in the 
first chapter (Oil and (fTBT) . 

If we have a bipartite state written as p = Y^nmuv Pm^\ m ) ® \f J )( n \ ® { v \-> 
where the Latin letters refer to the system A while the Greek letters refer to 
the system B, we can define the partial transposition operation (PT) p — ► p TB 
as the map which associate each density operator represented by p^ to the 
operator represented by p'mv So a partial transposition can be seen as a 
transposition made only with respect to the system B. 

Even if p is a physical state, the operator p TB could be unphysical (not 
positive semidefinite). This is a consequence of the fact that the transposition 
operation is a positive but not completely positive map. 

If p is a separable state, we have that the partial transposed is given by 



but since the operators p Bk correspond to physical states, than also the state 
of the whole composite system p is physical. 

So, observing that the trace and the Hermiticity of an operator are invari- 
ant under partial transposition, we can give the following criterion |15j . 

Theorem 8 (Peres-Horodecki). If a bipartite state p is separable, then the 
operator p TB is positive semidefinite. Conversely if p TB ^ 0, then the state is 
entangled. 

The condition p TB ^ is sufficient but not necessary for entanglement. In 
fact entangled states exist such that p TB > 0. These are the bound entangled 




(3.3) 



A: 
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states and cannot be distilled, while it can be shown that every state such 
that p TB ^ is distillable. 

Peres-Horodecki criterion is very good for discrete variables systems, like 
for example N qubits, while for continuous variables systems the following 
criteria are more suitable. 



3.2.2 Simon criterion 

Simon criterion is a CV adaptation of the Peres-Horodecki criterion. The basic 
idea is to consider the PT operation as a time reversal of a subsystem, so that 
in the phase space it is equivalent to a momentum inversion of system B. Let 
us give a proof of this. It can be shown that the Wigner function of a bipartite 
state can be written also in this form 

W(0 = X 5 J(q- x\p\q + x)e- i2xT Pd 2 x (3.4) 

where £ = (qA,PA,qB,PB) T , q = (<2X<?b) T and p = (pa,Pb) T - If we consider 
the representation of the density operator in the position basis, for the defini- 
tion of the partial transposition operation, we have that (qA, qB\p TB \q'Ai Qb) = 
(qa, QbIpWai 1b)- ^ e t us wr ite the Wigner function of p Ts : 

W(0 — ► W(0 = ^ J(q- x\p TB \q + x)e- i2xTp d 2 x. (3.5) 



Now, if we make the change of variables x — ► Zx in (|3.5p . where Z = 
diag[l, —1], we obtain 

W'(0 = J(q- x\ P \q + x)e~ l2xTZp d 2 x. (3.6) 

So we can conclude that the partial transposition acts as an inversion of the 
momentum ps and, if we define the PT matrix A = [1, 1, 1, —1], we can write 

PT : W(f) — ► W(A£). (3.7) 

The effect on the correlation matrix of a Gaussian state is then given by a 
matrix multiplication 

PT : V — ► AVA. (3.8) 

Using the result given in 112. 22j) . we can finally apply the Peres-Horodecki 
criterion to a CV system. 

Theorem 9 (Simon). If a bipartite state is separable, then 

AVA + in>0. (3.9) 
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Conversely if NVk + iVL has some negative eigenvalues, then the state is 
entangled. 

Moreover the non positive PT condition is also necessary for bipartite 
entanglement of Gaussian states. 

3.2.3 Duan et al. criterion 

Another criterion, which is a good test for CV entanglement has been proposed 
by Duan, Giedke, Cirac and Zoller [T7]. First of all we define two EPR like 
operators, which depend on the real parameter a, 

L a = \a\XA + -xb, M a = \a\pA - ~Pb, a£l. (3.10) 
a a 

Theorem 10 (Duan et al). If a bipartite state is separable, then 

VaGM, a/0 AL 2 + AM 2 > a 2 + 4r- (3-11) 

a z 

Dimostrazione. The inequality (|3.1ip can be derived supposing p in the sepa- 
rable form given in (|3.ip and then applying the Cauchy inequality. □ 

Another similar but weaker necessary condition for separability has been 
given by Mancini, Giovannetti, Vitali and Tombesi \18\ I19j. 

Theorem 11 (Mancini et al.). If a bipartite state is separable, then 

Vo€l, a / AL 2 a AM 2 > ^(a 2 + -^) 2 . (3.12) 

If we consider bipartite Gaussian states, the Duan et al. criterion can 
rearranged in order to give a necessary but also sufficient condition for separa- 
bility. In fact they shown that the CM of a bipartite state can be transformed 
by local symplectic operations to what they called standard form II, 



Vn 



( n\ —d\ \ 

n 2 d 2 
—d\ m\ 
\ d 2 m 2 ) 

where all the coefficients are positive and satisfy the following constraint 



(3.13) 



(mi - l)/(m - 1) = (m 2 - l)/(n 2 - 1) = {d 1 - d 2 ) 2 /{ ni - n 2 ) 2 = A 2 , A G R. 

(3.14) 
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Theorem 12 (Duan et al. - Gaussian states). A bipartite Gaussian state is 
separable if, and only if, when it is transformed to the standard form II we 
have 

AL\ + AMI > A 2 + 4, (3.15) 
where A is the parameter defined in j3.14\ )- 

3.3 Entanglement measures 

In the previous section we have seen several separability criteria which answer 
to the question: "Is this state entangled?". However, in the Quantum Infor- 
mation field, entanglement is a precious physical resource which one needs to 
quantify, like energy or entropy. So another question arises: "How much is 
this state entangled?". To quantify the entanglement of a state we need an 
entanglement monotone or entanglement measure, which is defined as a func- 
tion p — ► E{p) G R + which associates to each density operator a real pos- 
itive number and satisfies (possibly all) the following natural and reasonable 
properties. 

1 (Monotonicity) . The average entanglement of all the outcomes of a local 
operation and classical communication (LOCC) performed on p must 

be less than the entanglement of p. If pi are the results of a LOCC 

LOCC 

performed on p expected with probabilities pi, that is p ► p = 

Y,iPiPh then 

E(p)>J2piE(Pi)- (3-16) 

i 

2 (Convexity). Entanglement can not be increased in average by mixing 
other entangled states, i.e. 

E(J2PiPi) < I>i£(Pi)- (3-17) 

i i 

3 (Unitary equivalence). Entanglement must be invariant under local 
unitary operations, 

E(Ui U 2 • • • <g> U N pu\ ®U\---®U ] N ) = E(p). (3.18) 

4 (Faithfulness). E(p) = iff p is separable. 



5 (Additivity). E(jp\ (g> p 2 ) = E{p{) + E(p 2 ). 
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If we put together the quite abstract conditions 1 and 2, we get a very 
concrete physical requirement for an entanglement measure: 

E{p) > E(e(p)), (3.19) 

where p — ► e(p) is a LOCC. This means that an entanglement measure can 
not increase under local operations and classical communications. 

We are going to introduce three of the most important entanglement mono- 
tones: the entropy of entanglement, the (log-)negativity, and the entanglement 
of formation. 

3.3.1 Entropy of entanglement 

Definition 12. The entropy of entanglement is a measure defined for pure 
states. Given a bipartite pure state p = IV'KV'I £ 'Ha'&'H-b, then the entropy of 
entanglement is given by the Von Neumann entropy of one subsystem (tracing 
out the other), 

E H {p) = H[tr A (p)} = -tr[ti A (p) ln(tr A (p))]. (3.20) 

This is a consequence of the general property of quantum mechanics, ac- 
cording to which the more a quantum state is entangled the less we know about 
the single subsystems, so that the amount of disorder of a single subsystem 
quantifies the entanglement of the global system. 

Every bipartite pure state admits a Schmidt decomposition like 

|^) = ^c i |a i )®|6 i ), (3.21) 

i 

where {a{\ and {b{\ are two orthonormal basis of the respective Hilbert spaces 
TLa and TLb, and Cj £ M + are the Schmidt coefficients satisfying c| = 
1. The squared numbers cf are the eigenvalues of the reduced states tr^(/5) 
and tre(/3), which have the same spectrum. This implies that the entropy 
of entanglement is equal if computed on the subsystem A or B and it is a 
function of the Schmidt coefficients only, in fact 

E H (p) = H[tr A (p)} = H[ti B (p)} =~J2 c i lnc l ( 3 - 22 ) 

i 

The entropy of entanglement satisfies all the previous properties of an entan- 
glement monotone only if the initial state is pure. Moreover, it can be shown 
that if we require also a (weak) continuity property, then the entropy of en- 
tanglement is the unique entanglement monotone for pure states, in the sense 
that any other entanglement measure is a monotonic function of Eff(p). 
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An operational meaning of Eh{p) is connected with the distillable entan- 
glement. Given N copies of a state p = then, in the asymptotic limit 
N — > oo, one can distill NEh(p) copies of maximally entangled states. 

3.3.2 Entanglement negativity and logarithmic negativity 

Definition 13. Given a bipartite state p G TLa ®Hb, then the entanglement 
negativity is given by 

M[p) = \\p Tb ^~\ (3.23) 

where || . . . \ \t r is the trace norm defined as ||-A||t r = tr(v 2tA). 

The trace norm is equal to the sum of the singular values of a matrix, 
but in our case the density operator is Hermitian and therefore the singular 
values are equal to the modulus of the eigenvalues of p. Using the fact that 
tr(/5 TB ) = 1, eq. (|3.23p can be written also as 

Af(p) = ^ M, (3.24) 

Ai<0 

where Aj are the eigenvalues of p Ta , and the sum is made only over the negative 
ones. 

Now it is clear why N{p) is connected with the entanglement of p, in fact 
we have already seen that the condition p TB > is necessary for separability 
and so the negativity measures how much this condition is violated. 

Another entanglement monotone, which is not convex but it is additive, is 
the logarithmic negativity defined as 

E N (p) = ln\\p TB \\ tr . (3.25) 

It is an upper bound to the distillation rate, but a clear operational meaning of 
the log-negativity is not known. However En(p) is frequently used in Quantum 
Information, because it is easily computable if we deal with qubits or CV 
Gaussian states. 

Theorem 13. If p is a bipartite Gaussian state of 1x1 modes, characterized 
by its correlation matrix V , then 

E N (p) = max{- My), 0}, (3.26) 

where u is the minimum symplectic eigenvalues of the PT matrix KVA. 
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Dimostrazione. The matrix AVA is positive definite and therefore, by the 
Williamson theorem, it can be transformed with unitary operations to the 
diagonal form diag{s\, s±, S2, S2) given in (12.290 . The corresponding density 
operator is a tensor product of two thermal-like states p = p m ® p n2 with mean 
excitation numbers rtj = (s, — l)/2. Since AVA > then m > —1/2, moreover 
if n-i > the state is a physical thermal state, therefore the only possibility to 
have negative eigenvalues is when —1/2 < m < 0. The Fock basis expansion 
of a thermal-like state is 

^ = P.27) 

and if we suppose —1/2 < rtj < 0, then the trace norm is by definition 
IIAiJItr = Sfel n i| fc /(1 ~ Kl) fc+1 = (1 + 2n i )~ 1 = sT" . The conservation 
of the determinant det V = det(A^A) imposes that only one PT symplectic 
eigenvalue v can be less then one, so that Ej^(p) = max{— ln(i/), 0}. □ 

The log-negativity like any other entanglement monotone is invariant un- 
der local unitary operations. This means that the minimum PT symplectic 
eigenvalue v of a bipartite Gaussian state depends only on the symplectic 
invariants of the matrix V. In fact, if we write V in a block form like 



V = „ , (3-28) 



it can be shown that 




- ^ (3.29) 



where S = det A + det S - 2 det C. 



3.3.3 Entanglement of formation 

Definition 14. Given a quantum state p, the entanglement of formation is 
defined as 

E F (p) = minJ2 PkEH^k)^]) (3.30) 
fe 

with the constraint 

P = ^2pk\ipk)(ipk\- (3-31) 
k 

The operational meaning of Ep{p) is the minimum averaged entropy of 
entanglement of pure states, which are required to create the state p. 
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Since the decomposition (|3.3ip is not unique, the entanglement of forma- 
tion is difficult to calculate. A simple result is known for symmetric bipartite 
Gaussian states (with det A = det-B), for which we have that 

E F (p) = EhWvMA) (3.32) 

where \ifi r ) is a two modes squeezed state with the same entanglement nega- 
tivity of p. That is \ip r ) = X^n* annn ( r )/ cosh(r)|n)|ro), with r = Ej^(p)/2. 

A general result valid also for non-symmetric Gaussian states has been 
recently proposed [54] . 



Capitolo 4 

CV Teleportation 



Quantum teleportation is the transfer of an unknown quantum state form a 
sender (Alice) to a receiver (Bob) by means of a shared bipartite entangled 
state and appropriate classical communication. If the shared entanglement is 
infinite than Bob recovers an exact copy of the state sent by Alice. 

The first quantum teleportation protocol was given by Bennett et al. in 
1993 [11] for discrete variables systems. In this chapter we will consider only 
the continuous variable version, firstly proposed by Vaidman (1994) [12] and 
Braunstein et al. (1998) [T5] . 

The word "teleportation" could be misleading, since there is not a transfer 
of matter, energy or information, but only the transfer of a quantum state. For 
this reason, quantum teleportation is not in contradiction with the relativistic 
principle according to which a signal can not travel with a velocity greater 
than the speed of light and it is consistent with the no-cloning theorem [14J, 
which forbids to create a duplicate of a quantum state. 

From the point of view of the foundations of quantum mechanics, telepor- 
tation is not surprising in itself since it is simply an application of the more 
general exceptional phenomenon which Einstein called spooky action at a dis- 
tance. In fact we have already seen (EPR/Schrodinger paradoxes) that, if we 
act on a single part (Alice) of an entangled state, then we can change the 
state of the other (Bob). Now, if Alice and Bob, using only LOCC, make a 
proper distant action which project Bob state into Alice input state, then this 
is called quantum teleportation. 

In this chapter we first introduce the Braunstein-Kimble protocol using 
an infinitely squeezed EPR state and successively we give compact analytical 
results for the fidelity (success of the teleportation), in the situation when 
the shared bipartite state is a general Gaussian state which is not infinitely 
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entangled. 

4.1 Braunstein-Kimble protocol with infinite entan- 
glement 

We give a schematic description of the Braunstein-Kimble teleportation pro- 
tocol in the Heisenberg picture. 

1. Initial conditions 

Alice and Bob share the two modes of an infinitely squeezed EPR pair, 
such that 

Xb + X a = Pb -Pa = 0. (4.1) 

Alice has the input state p\ n described by the quadratures Xi n e pi n , that 
she wants to teleport. 

2. Beam splitter 

Alice mixes her mode of the EPR pair and the input state, through a 
symmetric beam splitter. The effect on the annihilation operators of the 
two modes is the following unitary transformation 



/ o + \ 1 / 1 1 \ / a in 
\ a- J ~ V2 \ 1 -1/1 a 

so the new quadratures become 



(4.2) 



(I 



a 3sin i %a - Pin ^ Pa / a n\ 

X± = - 7 =-, P± = — pr . (4.3) 

3. Bell measurement 

Alice makes two homodyne measurements of the quadratures x + and j5_ . 
If we call the results of the measured quadratures with two real numbers 
x+ and p— , then the system after the measurement is such that 

x a = -x in + y/2x+, p a = +p in - V2p— (4.4) 

Due to the entanglement correlations expressed in (|4.ip . the measure- 
ments made by Alice cause the following instantaneous collapse of the 
Bob mode 

Xb = X in ~ V2x + , p b = Pin ~ V2p-- (4.5) 
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4. Classical communication Alice sends to Bob, with a classical channel, 
the results of her measurements: x+ and p_. Without this classical 
communication Bob has a completely indefinite knowledge of the state 
of his mode. 

5. Conditional displacement 

Bob, according to the values x+ and p_ given by Alice, makes on his 
mode a correction displacement: 

%b * %out — %b \/2X-|_ , r a n\ 

i- ( 4 - 6 ) 
Pb > Pout =Pb + V 2p_ . 

This completes the teleportation, since now the output mode is described 
by the same quadratures of the input state 

x out = %in > ,^ ^ 

Pout — Pin ■ 

The protocol is very simple, in fact a great advantage of using CV variables 
rather than discrete variables is that the teleportation can be realized using 
only two simple optical elements: a beam splitter and a homodyne detector. 



4.2 Real teleportation with Gaussian states 

If the shared entanglement is not infinite, then the output state of the tele- 
portation will be not completely equal to the input. 

4.2.1 BK protocol in the Heisenberg picture 

Now we are going to see the relation between the input and the output state 
of the teleportation when the shared state is not an ideal EPR pair, but a 
Gaussian state p with zero mean displacement. We suppose the input state 
pin to be Gaussian and so, since the operations of the protocol are GCP maps, 
also the output state p ou t will be Gaussian. 

With the same notation that we used for Gaussian Wigner functions (|2.2ip , 
the degrees of freedom are then the first moments di n , the correlations Vi n of 
the input state pi n and the CM V of p. While, what we want to calculate are 
the first and the second moments of p ut, that we call d ou t and V ou t- 

A general treatment of the teleportation (good also for non Gaussian 
states) can be given using the Wigner function approach [29] . However, since 
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we suppose to deal with Gaussian states, we will calculate d ou t and V ou t in the 
Heisenberg picture, where it is easier to see the physical meaning underlying 
the equations. 

Let us define three column vectors containing the quadrature operators of 
the three modes (input, Alice and Bob) involved in the teleportation protocol, 



u a = Ha , u b = . (4.8) 






We define the EPR vector that in the ideal situation is zero but in general is 
a pair of operators given by 

ue = u b + Zu a (4.9) 

where Z = diag(l,—l). The correlations between the elements of this new 
vector can be expressed in a 2 x 2 matrix = (AuEiAuEj + AuEjAuEi) 
and if we write V in block form 



y= L n h ( 4 - 10 ) 




then, from the definition (|4.9j) . it is straightforward to prove that 

N = ZAZ + CZ + C T Z + B. (4.11) 

We define also an analogous vector, which contains the quadratures measured 
by Alice, 

U m = Uin + Zu a . (4.12) 

If we subtract (|4.12p to (|Q]> . we have 

Ub = u in + u E - Um, (4.13) 

which is actually an identity, but it gives us a useful way of expressing the 
quadratures of the mode possessed by Bob in terms of the fundamental oper- 
ators which are involved in the teleportation. Up to now, we have not started 
the protocol yet. When Alice makes the two homodyne measurements, she 
actually measures the vector ue with infinite precision so that it will collapse 
into a pair of real numbers 

measurement _ ~,o r i i t\ 

u E >u E eR. (4.14) 

Then Alice, using a classical channel, sends to Bob the vector ue, and Bob 
performs the conditional displacement r = ue to his mode, so that at the end 
of the protocol, (|4.13p becomes 

protocol A „ „ r a tr\ 

Ub > U out = U in + U E - (4-15) 
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Thanks to the simplicity of the Heisenberg picture, we can see that the effect 
of the Braunstein-Kimble protocol is to cancel the vector u m in (|4.13p . while 
the vector ue depends on the EPR correlations shared by Alice and Bob. In 
the ideal situation (14. 1|) ue is exactly 0, while in general it will introduce some 
additional noise that will affect the fidelity of the teleportation. 

Since we supposed the bipartite state p to be centered in phase space, then 
(ue) = and so, from (14.150 we always have di n = d ou t- While, for what 
concerns the second moments, from (14.150 we have that 

V out = V in + N. (4.16) 

Therefore, the difference between the input and the output state is the positive 
definite matrix N, which should be as much as possible near to zero in order 
to have a good teleportation. 

4.2.2 Teleportation fidelity 

A measure of how much a quantum state is similar to another one is the 
fidelity. 

Definition 15. Given two states pi and fa, the fidelity is defined as 

1 2 



(4.17) 



tr V VhhVh 

and if one operator is pure, e.g. p\ = \ip)(ip\, then (14.170 reduces to 

nw)u>\,P2) = w 2 h/>>- (4.i8) 

The fidelity !F{pi,P2) is a number between and 1 and can be seen as 
the probability that measuring the system fa it collapses to the state p%. In 
particular, if both states are pure, p\ = \ip)(ip\, p2 = \4>){4>\, then the fidelity 
is equal to the transition probability between the two states |(^|</>)| 2 - 

A quantitative measure of the success of the teleportation is then given 
by the fidelity between the input and the output state and if we suppose the 
input state to be pure then we have 

? = {lpin\Pout\tpin) = tr(pinPout)- (4.19) 

Using the Parseval theorem (12. 130 . the trace in the last part of (14.190 can be 
written as a scalar product of characteristic functions 

J7= hj X™(v)Xout(v)dv. (4.20) 
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The characteristic functions have the Gaussian form given in (|2.2U|) . so that 
(|4.2Qj) is a Gaussian integral which can be solved analytically 

T = — f e -h TnT (y™+Vout)a nd = 2 (-4 2i) 

2 W V^lVin + Vout] 

If we substitute the expression of V ou t (|4.16p , we obtain a formula which is true 
for the genera^] teleportation of an input pure state with CM Vi n , by using a 
bipartite state characterized by the EPR noise matrix iV given in (14, lip . 

T = , 2 (4.22) 
v/det[2V- n + iV] 

If we want a measure of the quality of the teleportation in itself which 
is independent from the input state, we can use the fidelity of entanglement 
swapping |27j . which is given by 

2 

Tswav = Vffiv (4 ' 23) 

Finally we give an important theorem, which underlines the importance of 
quantum entanglement to realize a good teleportation. 

Theorem 14. The maximum fidelity achievable in the teleportation of a co- 
herent state, based only on a classical strategy (without entanglement), is 
T= 1/2, 



4.2.3 A simple example 

As a simple example of quantum teleportation we consider a coherent input 
state (Vin = I) and a shared two-mode squeezed state characterized by the 
CM 

v= ( /cosh ^ -Zsinhr \ (4M) 
\ — Zsinhr Jcoshr I 

From (|4.1ip . the noise matrix is equal to iV = 2e~ r I, and using (|4.22p we get 
the simple result 

T=— — . (4.25) 

1 + e~ r y ' 



1 li the state shared by Alice and Bob has zero displacement and Bob always performs the 
"correct" conditional displacement, the teleportation protocol is invariant under displace- 
ment transformation, i.e., all states with the same covariance matrix but different coherent 
components are teleported with the same fidelity [271 129] . 
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It is interesting to note that if we use (|3.29p to calculate the minimum sym- 
plectic eigenvalue of the PT state, we get v = e~ r and so, we can write 

JT=— !— . (4.26) 

This is a particular example of the relation between the fidelity and the 
entanglement negativity which will be studied in detail in the next chapter. 

Finally we observe that the expression (|4.26p is consistent with (Thm. 
I14p . in fact if T > 1/2, then v < 1 and this is a sufficient condition for the 
entanglement of the two-mode squeezed state. 
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Capitolo 5 



Optimal fidelity and 
entanglement m 



5.1 Introduction 

The problem afforded in this chapter is: "Given the entangled state p shared 
by Alice and Bob, if we are allowed to perform local operations and classical 
communications, which is the maximum of the teleportation fidelity we can 
get?" and "What is its relation with the entanglement of pT\ The problem 
is non-trivial because the set of operations that Alice and Bob can adopt is 
very large. In fact, apart from local TGCP maps, they can adopt two further 
options: i) use non-trace preserving Gaussian operations in which some ancil- 
lary mode is subject to Gaussian measurement, i.e., projected onto a Gaussian 
state, rather than discarded [34J; ii) use local non-Gaussian operations (either 
with measurement on ancillas or not), i.e., those involving interactions which 
are non-quadratic in the canonical coordinates. The first class of maps, to- 
gether with TGCP maps, forms the most general class of Gaussian completely 
positive (GCP) operations, capable of preserving the Gaussian nature of the 
state shared by Alice and Bob. Non-Gaussian operations instead will trans- 
form the initial Gaussian bipartite state of Alice and Bob into a non-Gaussian 
one, and they can also increase the fidelity of teleportation in some cases |35j . 
Here we restrict to input coherent states, which represent the basic resource 
for many quantum communication schemes and because the conventional tele- 
portation protocol of Ref. |13j . while working excellently for coherent states, 
is less suited for teleporting nonclassical states [24] . 

The improvement of the teleportation of coherent states by means of local 
operations and its relation with the entanglement of the shared bipartite state 
has been already discussed in a number of papers [25j ESI E3 GS] . Ref. [25] 
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showed that in some cases the fidelity of teleportation may be improved by 
local squeezing transformations, while Ref. [26] showed that in the case of 
a shared asymmetric mixed entangled resource, teleportation fidelity can be 
improved even by a local noisy operation. These results were generalized in 
Ref. |27) which showed how fidelity can be maximized over all local trace- 
preserving Gaussian completely positive (TGCP) maps, i.e., those that can 
be performed by first adding ancillary systems in Gaussian states, then per- 
forming unitary Gaussian transformations on the whole system, and finally 
discarding the ancillas. Ref. [27\ confirmed that the optimal local TGCP map 
maybe a noisy one, i.e., that teleportation fidelity can be increased even by 
decreasing the entanglement and increasing the noise of the shared entangled 
state. Ref. [27], however, did not discuss the relationship between entangle- 
ment and the optimal fidelity J- op t- Ref- [23] instead found this relationship, 
but only for a subclass of symmetric Gaussian entangled state shared by Alice 
and Bob: for this class it is J- op t = (1 + > where v is the lowest symplectic 
eigenvalue of the partial transposed (PT) state (|3.29p .which is connected with 
the entanglement log- negativity by Ej^ = max[0, — lni/]. 

In this chapter we generalize in various directions the results of Refs. [27\ 
I23j . We show that if Alice and Bob share a bipartite Gaussian state with a 
given v and one restricts to local GCP maps which preserve such a Gaussian 
nature, the optimized fidelity always satisfies 

1+ " <F op t<^— ■ (5.1) 



1 + 3i/ - - 1 + v 
We also show that the upper bound is reached iff Alice and Bob share a 
symmetric entangled state. Moreover we determine the optimal local trans- 
formations at Alice and Bob sites and the corresponding value of J- op t as a 
function of the symplectic invariants of the shared CV entangled state when 
one restricts to local TGCP maps. 

The starting point of our investigation is the formula for the fidelity that 
we derived in the previous chapter 

T= - 2 (5.2) 
^det (2V in + N) 

where 

N = ZAZ + ZC + C T Z + B, (5.3) 

with Z = diag(l, —1). As we have already said, we shall restrict to the case 
of input coherent states, Vi n = I, so that Eq. (15. 2D reduces to 

T= -. 2 (5.4) 

v / 4 + 2TriV + detiV 
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The maximization of the teleportation fidelity over all possible Gaussian LOCC 
strategies therefore means to determine the optimal local transformation of 
matrices A, B and C which makes ./V as small as possible. 

As showed in [30], using an unbalanced beam splitter is equivalent, for 
the teleportation protocol, to a squeezing operation by Alice. Therefore the 
optimization over all Alice and Bob local operations includes also any eventual 
modification of the beam splitter used for the joint homodyne measurement. 

5.2 Upper and lower bounds for the fidelity of tele- 
portation 

In this section we prove Eq. (15. If) . i.e., the upper and lower bounds for the 
fidelity of teleportation for input coherent states. An important preliminary 
result enabling us to derive the two bounds is the fact that the optimal noise 
matrix N is very simple: in fact, the maximum teleportation fidelity is ob- 
tained when N is proportional to the 2x2 identity matrix I. More precisely, 
we have the following 

Lemma 1 (Optimal noise matrix). If u op t is an optimal local GCP map 
which gives the maximum of the fidelity J- max for the teleportation of a coherent 
state, then the resulting noise matrix is a multiple of the identity, that is 
N op t = 2n op tI. 

Dimostrazione. First of all we observe that for a 2 x 2, symmetric and pos- 
itive semidefinite matrix like N, the condition N = 2n op tI is equivalent to 
TriV = 2Vdet N. Therefore we have to show that oj op t is such that TrN op t = 
2y / det N opt . We do this by reductio ad absurdum supposing that u) op t gives 
a noise matrix with Tr N op t > 2-y/det N op t. However, within the class of local 
GCP maps, there exists a subclass of local symplectic (i.e., unitary Gaussian) 
maps realized by a generic symplectic Sb on Bob mode, and the associated 
symplectic map S a = ZSbZ on Alice mode, which act as an effective symplec- 
tic transformation on N opt , N' t = SbNoptS^ (see Eq. (15. 3p ). We can always 
choose Sb such that N' opt = -y/det N opt I, for which TrN' t = 2y / det N opt < 
TvN opt while det N' t = det N opt . However, we see from Eq. (|5.4|) . that this lo- 
cal symplectic operation increases the teleportation fidelity, but this is absurd 
because we assumed from the beginning that N op t is optimal. □ 

From this lemma and Eq. (15.41) we can therefore rewrite the optimal fidelity 
of teleportation in terms of the single positive parameter n op t = 
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as 

T< opt = — • (5.5) 

1 + n opt 

We can now derive an upper bound for T op t for a given entanglement of the 
state shared by Alice and Bob. We quantify such entanglement in terms of 
the lowest partially transposed (PT) symplectic eigenvalue v (I3.29p . 

Theorem 15 (Upper bound). For a given Gaussian bipartite state shared 
by Alice and Bob, with lowest PT symplectic eigenvalue v, the fidelity of the 
teleportation of a coherent state is limited from above by 

Topt < (5.6) 

Dimostrazione. Let us suppose that we can achieve a larger fidelity T = 1/(1+ 
n opt) with < n op t < v. Alice can in principle have at her disposal a two- mode 
squeezed state, with the usual correlation matrix 

(Jcoshr — Z smhr \ 
(5-7) 
— Zsinhr Jcoshr / 

(r is the squeezing parameter) and use this two-mode squeezed state, together 
with the bipartite state shared with Bob already optimized over all local GCP 
maps, to implement a CV entanglement swapping protocol [51] . In fact, by 
mixing at a balanced beam splitter her mode of the bipartite state shared with 
Bob and one part of the two-mode squeezed state, and performing homodyne 
measurements at the output, Bob mode gets entangled with the remaining 
part of the two-mode squeezed state in Alice hands. Since the noise added to 
the teleported state is N op t = 2n op tI, it is straightforward to see that the two 
remaining modes are then described by the following CM 

(/coshr — Zsinhr \ 

y ■ 1 m , , , ■ (5-8) 

— Z sinn r 1 [2n op t + cosh r\ J 

In other words, before entanglement swapping, Alice and Bob shared an en- 
tangled state with CM V and entanglement characterized by v; after entan- 
glement swapping, they share a state with CM W swap . In the limit of infinite 
squeezing the lowest PT symplectic eigenvalue of W swap tends to n opt , i.e., 
lim r ^oo inswap = n opt . Since we supposed n op t < v, this means that for a 
sufficiently large squeezing parameter r, ^ su , ap < u, i.e., Alice and Bob have 
increased their entanglement. However this is impossible because we have 
employed only local operations. Therefore it must be n op t > v. □ 



Chapter 5 - Optimal fidelity and entanglement 



37 



We complete the characterization of the optimal fidelity of teleportation 
in terms of the entanglement shared by the two distant parties by providing 
also a lower bound for J- op t, proving in this way the result of Eq. (15. 1|) . 

Theorem 16 (Lower bound). For a given Gaussian bipartite state shared 
by Alice and Bob, with lowest PT symplectic eigenvalue v, the fidelity of the 
teleportation of a coherent state is limited from below by 

Dimostrazione. From the definition of symplectic eigenvalue, one has that a 
4x4 symplectic matrix S exists which diagonalizes AVA (A = diag(Z, I)), i.e., 
the PT matrix of the CM V. This means SAVAS T = diag(^, u, where 
\i is the largest PT symplectic eigenvalue. By writing S in 2 x 2 block form 

S-( W ' W >), (5.10) 

and rewriting the diagonalization condition for the upper 2x2 block only, one 
gets the following condition 

W a ZAZWj + W a ZCW^ + W h C T ZWl + WbBW? = vl. (5.11) 

The symplectic transformation S transforms the vector of quadratures £ into 
I' = (x' a ,y' a ,x' b ,y' b ) T = Si and the PT vector into |» = (x^,x£,y£) T = 
SA£. One has [x' a ,y^] = i, because commutation relation are preserved by S, 
implying 

det W a + det W b = 1. (5.12) 

The commutation relation is instead not preserved for the PT transformed 
quadratures, and introducing a real parameter e such that = ie, we 

get another condition for the two upper blocks of 5*, 

- det W a + det W b = e, (5.13) 

which together with Eq. (|5.12p . gives the parametrization 

det W a = (l-e)/2, det W b = (1 + e)/2. (5.14) 

Now, since Aa^ 2 = Ay" 2 = u/2, the Heisenberg uncertainty principle imposes 
that |e| < v and in particular for every entangled state we have |e| < v < 
1. This latter condition, together with Eq. (15.140 . suggests an alternative 
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parametrization in terms of the angle 9 = arctan v/(l — + e) (0 < 9 < 
tt/2), 

VdetWa = sin 9, y/detW b = cos9. (5.15) 

The 2x2 matrices W a and W b and the parameter 9 allow to construct an 
appropriate local map which will lead us to derive a lower bound for the 
fidelity. This local map is a TGCP map which, at the level of CM, acts as 



V = SVS T + G, (5.16) 

with S and G satisfying 

G + ij - iSJS T > 0. (5.17) 

If the TGCP map is local, then S = S a © S b and G = G a © G b , with G k + 
iJ — iSkJSfr > (k = a,b). 

The desired local TGCP map ujg is defined in terms of S a , S b , G a and G b 
in the following way 

_ ( WS-r 0< e <,/4 
| ZWaZlsin^]- 1 vr/4<#<^/2, 

* = / °<^"/ 4 (5 , 8b) 

\ Wbfsin^]" 1 vr/4 < < tt/2, 

= J [1-tan^]/ 0<^<vr/4 

\ tt/A<9< tt/2, K ' 7 

f 0<6 , <tt/4 , k 

G b = < r . n ~ 7 5.18d) 

\ [l-cot 2 0]/ tt/4<^<vt/2- 

By applying Eqs. (|5.3|) . (|5.11|) and (|5.16|) . one can see that this local TGCP 
map transforms the noise matrix N into a final matrix proportional to the 
identity, given by 

N = [v/ cos 2 6 + 1 -tan 2 9]I, < 9 < it/A, (5.19) 
N = [u/ sin 2 9 + 1- cot 2 9}I, tt/A<9<tt/2. (5.20) 

It is however convenient to come back to the parametrization in terms of e, 
which allows to express the final N in a unique way, for < 9 < tt/2. In fact, 
from Eqs. ([CTJ]) - (|520J ). one gets 



(5.21) 
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which, inserted into Eq. (|5.5p . yields 



i + kl 



l + i/ + 2e 



(5.22) 



From the condition imposed by the Heisenberg uncertainty principle < |e| < 
u, we see that the fidelity is minimum when |e| = v, so that we get the following 
lower bound 



Ft opt > 



1 



1 + 3^ 



(5.23) 
□ 



Theorems 1 and 2 provide a very useful characterization of the optimal 
fidelity which can be achieved with Gaussian local operations at Alice and 
Bob site. In fact, the bounds are quite tight because the region between the 
upper and the lower bound is quite small (see Fig. 15. Therefore, by simply 
computing the lowest PT symplectic eigenvalue of the CM of the shared state 
and using the bounds, one gets a good estimate of the maximum fidelity that 
can be obtained with appropriate local operations. In fact, the error provided 
by the bounds is never larger than 0.086 (see Fig. 




Figura 5.1: Plot of the upper and lower bounds (Eq. (|5.6p and (|5 .9 j) respec- 
tively) for the fidelity of teleportation of coherent states. The blue region is 
the allowed region in the (J-, v) plane. 



Corollary 1 (Upper bound achieved in the symmetric case). The upper 
bound J- opt = 1 /(l + l/ ) *s achieved iff the bipartite Gaussian state shared by Al- 
ice and Bob is symmetric. The optimal local transformation in the symmetric 
case is a local symplectic map. 
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Figura 5.2: Plot of the distance between the upper and lower bounds for 
the teleportation fidelity versus the allowed values of v for an entangled state 
between Alice and Bob. We see that the error AJ- under which we can estimate 
the maximum fidelity is less than 0.086. 



Dimostrazione. The "if" part of the theorem directly follows as a special case 
of the preceding proof. If the Gaussian state shared by Alice and Bob is 
symmetric, it is det W a = det W b , implying e = 0. Then, Eq. (|5.22|) shows that 
in this case the fidelity reaches the upper bound, J- op t = 1/(1 + v). Moreover 
in this case = 7r/4 and the local TGCP map of Eqs. (|5.18p is optimal and 
it is a symplectic one, with S a = ZW a zV2, S b = W b V2, G a = G b = 0. The 
"only if" part instead can be easily proved by using the result of Theorem 3 
about the CM of the optimized bipartite state shown in the following section. 
The proof is given in the Appendix. □ 

This latter corollary provides the generalization of the result of Ref. [23J , 
which obtained the same relation between optimal fidelity and v but by consid- 
ering only a special class of symmetric Gaussian bipartite state for Alice and 
Bob, obtained by mixing at a beam splitter two single-mode thermal squeezed 
states. 



5.3 Determination of the optimal local map 

We have derived a lower bound for the optimal fidelity of teleportation of 
coherent states, by explicitly constructing the family of local TGCP maps 
ujq of Eq. (j5.18j) . which transform Alice and Bob shared state so that the 
corresponding fidelity of teleportation is given by Eq. (|5.22p . interpolating 
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between the lower and upper bound of Theorems 1 and 2 by varying e = cos 29. 
The map ujq is symplectic only for 9 = tt/4 and in this case it is the local 
optimal map for a symmetric shared state, since it reaches the upper bound 
of Theorem 1 (see Corollary 1). When 9 ^ vr/4, ujq is a noisy (i.e., non 
unitary) map, and it is not optimal in general, because one cannot exclude that 
different LOCC strategies by Alice and Bob may yield a better noise matrix 
N and therefore a larger value of the teleportation fidelity. As discussed in 
the introduction, here we shall study the optimization of the teleportation by 
restricting to GCP maps, which preserves the Gaussian nature of the bipartite 
state initially shared by Alice and Bob. 

In this section we shall derive two results: i) the general form of the final 
CM of the bipartite Gaussian state after the optimization over all local GCP 
maps; ii) the optimal local TGCP map, i.e., the local TGCP map which 
maximizes the teleportation fidelity when one restricts to TGCP maps only, 
excluding in this way measurements of ancillary modes. 

Ref. [27] already provided the analytical procedure for the determination 
of all the parameters of the optimal TGCP map. Here, by further elaborating 
the approach of Ref. [27J, we will show that the optimal TGCP map can always 
be written in a simple form, as a local symplectic operation eventually followed 
by a single mode attenuation [36], either at Alice or Bob site. 

5.3.1 Standard form of the optimal correlation matrix 

In this subsection we show that, even if we do not know the specific form of the 
optimal GCP map, one can always characterize it indirectly by determining 
the general form of its outcome, i.e., the general form of the CM of the final 
Gaussian state shared by Alice and Bob after the maximization. We begin 
with the following lemma. 

Lemma 2 (Standard form III). The correlation matrix V of every bipar- 
tite Gaussian state can be transformed by local symplectic operations into the 
following normal form 

( n\ —d\ ' 

n 2 d 2 
—d\ mi 
\ d 2 m 2 J 

where all the coefficients are positive and satisfy the following constraints 

n\ — n 2 = mi — m i = d\ — d 2 = A, A£l. (5.25) 



(5.24) 
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That is: 



( n + X 



n 



-d — X 



-d-X 
m + A 



d 



\ 



(5.26) 



/ 



d m 

Dimostrazione. It is well known that it is possible to transform every V of an 
entangled state in the usual normal form (standard form I) |17t [T6] 

/ a —ca \ 



Vn 



-ci 



V 



-ci 

a c 2 

6 



C2 



(5.27) 



where all the coefficients are positive. Now we perform a local symplectic 
operation composed by two local squeezing operations, S a = diag(y / r^, 
and Sb = diag(- v /rfc, l/y^). We impose the first two conditions n\ — 112 = 
mi — m,2 = A, 

A, 



a(r a ~r a 1 )= b(r b - r b x ) 
which solved for positive r a and r& give 

r a (A) = A/(2a) + ^/l + (A/2a) 2 , 



r b (X) = X/(2b) + y/1 + (A/26) 2 . 
Now we impose the last constraint d\ — c?2 = A, that is 

ci \/ r a(A)r b (A) - c 2 / \/r a (X)r b {X) = X. 



(5.28) 

(5.29) 
(5.30) 

(5.31) 



Our lemma is proved if there is at least one solution A of Eq. (|5.3ip . Since 
A = is the trivial solution when c\ = C2 and Vn = V%, we can exclude this 
particular case and divide Eq. (|5.3ip by A. Therefore we have to show that 
the equation 

1 



/(A) = -Mr a {X)r b {X) - c 2 /yjr a {\)r h {\)] = 1 



A 



(5.32) 



admits at least one real solution. If |A| S> a and |A| ^> b, we can power expand 
the square roots in Eqs. (I5.29|) - ()5.30p so that we easily find the following limits 
for /(A): 



hm /(A) 

A— >oo 

hm /(A) 

A — > — OO 

lim /(A) 



ci/Vab < 1, 
c 2 /Vab < 1, 
±sign(ci — c 2 )oo. 



(5.33) 
(5.34) 
(5.35) 
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The inequalities in Eqs. (|5.33p - (|5.34|) follow from the Cauchy- Schwartz in- 
equality (a?!) (a;!) > {X1X2) 2 applied to the quadrature operators of the two 
modes. Given the three limits (|5.33j) . (I5.34|) and (|5.35j) . since /(A) is continu- 
ous everywhere except at the origin, at least one solution of Eq. (|5.32p exists. 
Moreover this solution has the same sign of c\ — c 2 - □ 

We have defined the standard form of lemma 2 as standard form III because 
it is very similar to the standard form II defined in Ref. [T7j for the determina- 
tion of a necessary and sufficient entanglement criterion for bipartite Gaussian 
states. In particular the two standard forms coincide in the special case of a 
symmetric bipartite state (n± = mi and n 2 = m 2 or equivalently n = m). 

Theorem 17 (Form of the CM of the optimized bipartite state). The 

optimal GCP map uj op t maximizing the teleportation fidelity is such that the 
CM of the transformed bipartite state is in the standard form III V\ defined 
by Eqs. fQ^-pjgp. 

Dimostrazione. By means of local symplectic operations, we can always put 
the CM of the bipartite state of Alice and Bob in the form of Eq. (I5.24p . but 
without the constraints of Eq. (|5.25p . We first restrict to local symplectic oper- 
ations and show that the optimal local symplectic operation always transforms 
to a state with a CM satisfying the constraints of Eq. (|5.25p . Since the CM 
is tridiagonal, then any possible optimal map must be a squeezing transfor- 
mations of the two modes given by S a = diag(r a ,r~ 1 ) and S b = diag(r;,, r^ 1 ) 
[2"7] . Let us define 



a{r a ,r b ) = r a ni - 2d 1 r a r b + rfmi, (5.36) 



P{r a ,rb) = r a 2 n 2 - 2d 2 (r a r b ) 1 +r b 2 m 2 , (5.37) 



so that the noise matrix of Eq. (15. 3D is equal to iV = diag(a, (3). The optimal 
map must minimize det iV = a/3, and therefore we impose that 

Va(r a ,r b )P(r a ,r b ) = 0, (5.38) 

where V = (d Ta ,d Tb ). Due to Lemma 1 , we must also have that ct(r a ,r{,) — 
0(r a ,rf,) 7^ 0, and therefore Eq. (|5.38p reduces to 

V[a(r a ,r b )+(3(r a ,r b )]=0. (5.39) 

The CM of the transformed state is the optimized one iff the optimal local 
symplectic operation is the identity map, that is, if a(l, 1) = (3(1, 1) and 

V[a(r a ,r b )+P(r a ,r b )} =0. (5.40) 

r a =r b =l 
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It is easy to check that these conditions are satisfied iff 

n\ — ni = mi — = d\ — d,2, (5-41) 

which are exactly the constraints of Eq. (I5,25p , The theorem is proved if we 
show that the normal form V\ of Eq. (|5.26p is actually kept also if one maxi- 
mizes over the broader class of GCP maps. In fact, if by reductio ad absurdum, 
we assume that an optimal (non-symplectic) GCP map exists leading to a CM 
not satisfying the constraints of Eq. (|5.25|) . we could always apply a further 
symplectic map which, by repeating the maximization above, would transform 
to a state with a CM satisfying the constraints (|5.25|) and yielding a larger 
teleportation fidelity. But this is impossible, because it contradicts the initial 
assumption of starting from the optimal bipartite state. □ 

In other words, the CM of the state shared by Alice and Bob after the 
maximization of the teleportation fidelity is the one with the standard form 
III of Eq. (|5.26p because it is the unique CM for which the optimal map is the 
identity operation on both Alice and Bob site. 

5.3.2 Optimal trace-preserving Gaussian CP map 

In the former subsection we have determined the form of the CM of the opti- 
mized state of Alice and Bob, without determining however which is the local 
GCP map which maximizes the teleportation fidelity. Here we find this optimal 
map, restricting however to the smaller class of trace-preserving GCP maps. 
The case of Gaussian maps including Gaussian measurements on ancillas will 
be afforded elsewhere. 

Ref. [32j has introduced the notion of minimal noise TCGP maps, as the 
extremal solution of the condition of Eq. (I5.17p . These maps are the ones 
that, for a given matrix S, possess the "smallest" positive matrix G realizing 
a CP map. It is easy to check that a minimal noise TGCP map satisfies the 
relation det G = (1 — det S) 2 . An example of minimal noise TGCP map is an 
attenuation [36j, i.e., the transmission of a single boson mode through a beam 
splitter with transmissivity r (0 < r < 1), such that 

a ra + ^fl-T 2 a v , (5.42) 

where a = (x + ip) jyf 7 ! is the annihilation operator of the mode, and ay that 
of the vacuum mode entering the unused port of the beam splitter. 

It is evident that the TGCP map maximizing the teleportation fidelity has 
to be a minimal noise TGCP map [27]. We prove now a useful decomposition 
theorem. 



Chapter 5 - Optimal fidelity and entanglement 



45 



Theorem 18 (Decomposition of TGCP maps). A minimal noise TGCP 
map uj on a single mode system with det G < 1 can always be decomposed into 
a symplectic transformation a\, followed by an attenuation r and by a second 
symplectic transformation gi, that is, 



Therefore a local minimal noise TGCP map on a bipartite CV system can al- 
ways be decomposed into a local symplectic map, followed by the tensor product 
of two local attenuations and by a second local symplectic map. 

Dimostrazione. We consider a generic minimal noise TGCP map such that 
V -»• V' = SVS T + G for a generic CM V, with det G = (1 - detS) 2 . G 
is a positive symmetric matrix, and therefore a symplectic matrix T2 exists 
such that G = T^T^i^ — s), where s = det S. We then define the symplectic 
matrix T\ = ■^T 2 ~ 1 5, and we also consider an attenuation map with trans- 
missivity yfs. If we now first apply the symplectic map defined by T%, then 
the attenuation map and finally the second symplectic map defined by T2, by 
using the relations T2y/sT\ = S and G = T2T 2 r (l — s), one can check that the 
composition of the three maps reproduces the given TGCP map. □ 

Corollary 2. A minimal noise TGCP map 10 on a single mode system with G 
proportional to the identity matrix, i.e., G = (1 — s)I (0 < s < 1) can always 
be decomposed into a symplectic transformation o\ , followed by an attenuation 
r, 



Therefore a local minimal noise TGCP map on a bipartite CV system with 
Gi = (1 — si)I (0 < Si < 1, i = a,b) can always be decomposed into a local 
symplectic map, followed by the tensor product of two local attenuations. 

Dimostrazione. It is sufficient to repeat the former proof and consider that, 
since G = (1 — s)I, T2 = I and therefore the second symplectic map is the 
identity operation. □ 

This latter case is of interest because the optimal TGCP map must have 
in fact the property Gi = (1 — i = a,b. To show this we first simplify 
the scenario by exploiting the results of Ref. [27], which provides the gen- 
eral analytical procedure to derive the optimal local TGCP map. In fact, 
Ref. [27] shows that when the optimal local TGCP map is a minimal noise, 
non-symplectic one, it can be performed on one site only, i.e., either on Al- 
ice or on Bob alone. Suppose that the non-symplectic map is performed on 



UJ = 0~2 T 0~\ ■ 



(5.43) 



uj = r o o\ . 



(5.44) 
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Bob site; it is straightforward to see that, under a generic local TGCP map 
V -> V = SVS T + G, with S = I ® S b and G = I © G b , the noise matrix N 
transforms according to N ^ N' = T + G b where 



r = ZAZ + ZCSl + SbC 1 Z + S b BS\ 



T r 



vr 

'6 • 



(5.45) 



Ref. [27] shows that the optimal map is such that r oc G b , and since Lemma 
1 shows that the optimal N is proportional to the identity, this implies that 
G b must be proportional to the identity. Therefore Corollary 2 leads us to 
conclude that the optimal local TGCP map is either a local symplectic map, 
or a local symplectic map followed by an attenuation by a beam splitter, placed 
either on Alice or on Bob mode. 

Theorem 4 and Corollary 2 therefore provide a very simple and clear de- 
scription of the TGCP map which maximizes the teleportation fidelity, which 
is not evident in the treatment of Ref. [27]. We can further characterize the 
optimal local TGCP map by determining: i) the form of the first local sym- 
plectic map; ii) the conditions under which the optimal local operation is noisy, 
i.e., when one has also to add a beam splitter with appropriate transmissivity 
on Alice or Bob mode in order to maximize the teleportation fidelity. In order 
to do that we first need a further lemma, similar to lemma 2. 

Lemma 3. For any given positive real parameter r/, the correlation matrix 
V of every bipartite Gaussian state can be transformed by local symplectic 
operations into the following normal form 



V„ 



( n\ —d\ \ 

n 2 d 2 
—d\ mi 
\ d 2 m 2 ) 



(5.46) 



where all the coefficients are positive and satisfy the following constraint 

n\ — n 2 = rj(di — d 2 ) = n 2 (m\ — m 2 ) = A, Agi. (5-47) 
That is, we have a family of normal forms depending on the parameter r\, 



( n + A 



Vn 



-d-X/rj \ 
n d 
d — A/17 m + -V 7 ? 2 

d m ) 



(5.48) 
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Dimostrazione. With the same procedure used in the proof of Lemma 2, we 
arrive to an equation similar to (|5.32p . while the corresponding three limits 
are exactly the same of (I5.33p . (|5.34p and (|5.35j) . since the factor 77 cancels 
out. As a consequence the continuity argument is valid also in this case, and 
therefore, for every fixed parameter 77, one can find a transformation which 
puts the CM in the normal form V n . □ 

We notice two facts that will be useful for the next theorem: i) the optimal 
CM standard form of Theorem 3, Vx, belongs to the class of normal forms V^, 
since it is obtained for 77 = 1; ii) when < 77 < 1, the state with CM 
is transformed into the Gaussian state with CM equal to V\ when a beam 
splitter with transmissivity 77 is put on Bob mode. We then arrive at the 
theorem about the optimal TGCP map. 

Theorem 19. The optimal local TGCP map maximizing the fidelity of tele- 
portation of coherent states can always be decomposed into a local symplectic 
map, eventually followed by an attenuation either on Alice or on Bob mode. 
The first local symplectic map is the one transforming the CM of the Gaussian 
state shared by Alice and Bob into one particular normal form of the family V v 
defined in Eqs. fr5.46\ )- [5-4'A ); with < 77 < 1. One has to add the attenuation 
on one of the two modes for realizing the optimal TGCP map if there is a value 
of r\, let us say r\ = r , such that the coefficients of V T satisfy the relations 

t = —^—, r<l. (5.49) 
m — 1 

// instead the condition of Eq. J5,^^| J is never satisfied in the interval 77 £ 
[0,1), the optimal TGCP map is formed only by the local symplectic map 
transforming the CM into the normal form V\ (i.e., 77 = 1 and no attenuation 
is required). 

Dimostrazione. Using Lemma 3, we can always apply a local symplectic map 
which transform the CM into one of the form of the family of Eq. (|5.48|) , with 
77 = r. We then apply an attenuation on Bob mode with transmissivity r and 
then try to find the maximum fidelity as in the proof of Theorem 3. We have 
that the two diagonal elements of the noise matrix N now read 

q (t) = /3(t) =n-2Td + T 2 m + l-T 2 , (5.50) 



The optimal map must minimize det N 

da(r) 
dr 



= a 2 , and therefore we impose 

= 0, (5.51) 
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which is satisfied iff r = d/{m — 1). If < r < 1, this map composed by the 
local symplectic map and the attenuation is the optimal map. If instead for 
any r G [0, 1], the condition of Eq. (15. 49ft is not satisfied, there is no critical 
point in this interval and therefore the optimal map is just the symplectic 
transformation to the normal form V\. One has also to check the behavior at 
the lower boundary value r = 0, but this is trivial because this means that Bob 
uses the vacuum to implement the teleportation which is never the optimal 
solution if we have an entangled channel. In fact, if Bob uses the vacuum, 
the channel looses its quantum nature and the maximum of the fidelity is the 
classical one T = 1/2, which is below the lower bound for any entangled state 
given by Eq. (|5.9p . □ 

Theorem 5 therefore characterizes in detail the optimal TGCP map, giving 
in particular the conditions under which this map is noisy, i.e., non-symplectic 
and therefore when teleportation is improved by increasing the noise and 
decreasing the entanglement of the shared state. 

Using this latter theorem we can also determine how, from an operational 
point of view, one can compute the value of the teleportation fidelity max- 
imized over all TGCP maps, starting from the symplectic invariants of the 
bipartite Gaussian state initially shared by Alice and Bob. From the CM of 
this latter state one can: 

1. compute the four symplectic invariants a = V det A, b = V det B, c = 
y 7 ! detC| and v = det V. 

2. Knowing the first three invariants of the channel, the elements n, m, and 
d of the normal form V v can be expressed as functions of only the two 
unknown parameters A and r\ as 



3. The two parameters A and r\ can be found solving the following system 
in the region < r\ < 1, 



n(A) 
m(A, n) 
d(\,rj) 



A/2 + y/a? + (A/2) 2 , 



(5.52) 
(5.53) 
(5.54) 



A/2r? 2 + y/b 2 + (A/27/ 2 ) 2 , 



A/27/ + Vc 2 + (A/2 ?? ) 2 . 




(5.55) 



and solving also the first equation in the boundary r\ = 1 



det Vi(A) = v, 



(5.56) 
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(the two conditions of (|5.55p come from the invariance property of the 
determinant of the channel and form the maximization condition of 
Eq. (Ml ). 

4. We call (Aj,r?i) with % = l,2...k, the union of the solutions of (15.551) 
and (15.560 (we have at least one solution because (15.560 admits at least 
a solution). We then compute the candidate fidelities 



so that the maximum fidelity will be J- op t = max{Fj}. 

5.4 Appendix 

We now prove the "only if" part of Corollary 1, i.e., that if T op t = 1/(1 + v) 
(the upper bound of the optimal fidelity), then the bipartite Gaussian state 
shared by Alice and Bob is symmetric. 

Dimostrazione. Theorem 3 shows that the final CM after any optimization 
map must be V\ of Eq. (|5.26j) . Using Lemma 1 and the explicit form of 
Vi, the hypothesis is equivalent to n + m — 2d = 2u, so we can make the 
substitution d = (n + m)/2 — u, which is just a different parametrization: 
Vx(n,m, d, A) — > Vi(n,m, u, A). Now, the condition that v is equal to the PT 
minimum symplectic eigenvalue gives us a constraint on the parameters of the 
matrix V\ 



If the state is symmetric, which means n = m, then the condition of Eq. (|5,58p 
is identically satisfied. If the state is non symmetric, Eq. (|5.58|) is a non- 
trivial equation that solved for A gives A = (m — n) 2 /8v — n — m. However 
the corresponding matrix V\ = (n, m, v , A) is not the CM of a physical state; 
in fact, the characteristic polynomial of V% can be written as P{x) = (cq + 
c\x + x 2 )(go + g\x + x 2 ) where cq = —v(n + m + v), but this means that V\ 
has at least one negative eigenvalue and therefore it is not positive definite. 
Therefore J- op t = 1/(1 + v) is realized only if Alice and Bob state is a Gaussian 
symmetric state. □ 



F t = 2 2 + ^a 2 + \ 2 /4+^b 2 V l + \ 2 /4 




(5.57) 



u{Vi(n,m, v, A)} 



(5.58) 
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Parte II 

OPTOMECHANICAL 
DEVICES 



Capitolo 6 

Quantum Optomechanics 



In this chapter we give some fundamental theoretical tools, useful to study the 
dynamics of optical and mechanical systems and their mutual interaction. 

There is one inconvenience that is usually neglected in abstract Quantum 
Information theories, but that we need to take into account when considering 
real systems; this is noise, thermal and quantum noise. Thermal noise, is the 
same which is studied in classical mechanics, is due to the thermal excitation of 
the system and vanishes when T — ► 0, while quantum noise is due to intrinsic 
quantum uncertainties of non-commuting operators and cannot be completely 
eliminated even at zero temperature. Both thermal and quantum noises can be 
described by different approaches: the master equation (Schrodinger picture), 
the Fokker-Planck equation (phase space), the Langevin equations (Heisenberg 
picture). Here we are going to consider only the latter approach, which is 
suitable for linearized systems (the ones we are interested in) and moreover it 
gives a consistent treatment of quantum Brownian motion. 

The last section is about the coupling between a mechanical resonator and 
an optical mode, due to the effect of radiation pressure. 



6.1 Input-output theory 

The dynamics of an optical cavity mode can be formulated using an input- 
output theory, in which the noise is seen as an input filed entering in the cavity 
or, equivalently, as an output field exiting from the cavity. 
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Figura 6.1: Optical cavity mode interacting with input and output fields. 



We derive the Langevin equations following the treatment of Gardiner and 
Collett [55]. We start form the full Hamiltonian 



H = H c + H, 



bath 



H 



hit ■ 



(6.1) 



where H c is the Hamiltonian of the cavity mode (we neglect constant vacuum 
terms) 

H c = hcOcO^a, (6-2) 

Hbath is the Hamiltonian of the external fields modeled as a bath of harmonic 
oscillators 



H 



bath 



duj hujtf (uj)b(uj) 



(6.3) 



while Hi n t is the linear interaction Hamiltonian in the rotating wave approxi- 
mation 



Hint = ih I duoK{uj) Mr(oj) — b{uo)a^ 
Jo L 



(6.4) 



where k(u) is the coupling constant between the modes. Consistently with 
the rotating wave approximation, we can assume that the integral in (|6.4p 
gives a time average non-zero contribution only for frequencies near the cavity 
resonance u> e ~ 10 15 , so that we can extend the lower integration limit to — oo, 



Hi, 



int 



/oo r 
cIlok,(uj) atf(uj) — b(ui)i 



The commutation relation are the usual bosonic CCR, 

[0,0+] = 1 (b{uj),b\J)]=8(u-J), 



(6.5) 



(6.6) 



from which it is evident that a and are dimensionless operators, while b(u>) 
and b'(oj) have the dimensions of a square root of time. Let us apply the 
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Heisenberg equation of motion, given in the first chapter (|1.12[) . to the bath 
annihilation operator, 

b(u) = —\b(u),H] = -iub(u) + k(uj)cl, (6.7) 
in 

which integrated for t > to gives 

S(c) = e'^-^boiu) + k(u) f dt'e-^-^ait'), (6.8) 

J to 

where bo(u>) is the bath operator at the initial time to. We can write the same 
solution for t < t\ 

b(u) = e-^-^bxiu) - k{lo) / dt'e-^- 1 } a(0, (6.9) 



where bi(u) is the bath operator at the final time t±. 

Now we consider the Heisenberg equation for the cavity annihilation oper- 
ator, 

1 f°° 
a = — [a, H] = — iuj c a — I dwK[uj)b(oj), (6.10) 

J — oo 

if we substitute the solution for b{uj) depending on the initial condition (16.8p . 
we get 



a = —iuj r a 



dcuK{uj)e- iuj{t - to) b {oj) - / dj/e 2 (w) / dt'e^-^ait'). 

) J — oo J to 

(6.11) 

Now we make the Markovian assumption of a frequency-independent coupling 
constant k(u) = 7/27T, so that (|6.1ip becomes 



due- iuj{t - to) b (uj) - 7 / dt'S(t-t')a(t'). (6.12) 
) J t 

We define the input-output field operators 

i r°° 

fiin(t) = -7T / dtue-^-^kitu), (6.13) 



2?r . 

i r°° 

aout(t) = — / due-^-^hiu), (6.14) 

satisfying the CCR 

[ain(t),ai(t')\ = [d out {t),d\ ut {t')] = 6(t - t'). (6.15) 

If we apply the property of the delta function according to which, if it is 
centered on one of the integration limits, it gives only half of its contribution 

b dx5(x-b)f(x) = \f(b), (6.16) 
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then (|6.12p becomes 



-iu c a - -a+ yfya in , (6.17) 



while if we repeat the same procedure using the final time solution for b(u) 
dSHD, we get 

a = -iu} c a + -a - y/^a out . (6.18) 
The difference between (|6.18p and (|6.17p give us a useful input-output relation 

hn + a out = y^a. (6.19) 

Equations (|6.17p and (|6.18p are two equivalent quantum Langevin equations 
describing the cavity mode dynamics with a damping rate 7 and where the 
noise field is seen as an input or output optical field. However, since a ou t 
depends on future unknown boundary conditions, Eq. (|6.18|) is useless form a 
practical point of view, while (|6.17p gives us the possibility to easily calculate 
the time evolution of the first and the second moments of a, thanks to the 
statistical properties of hi n . 

In fact, if at the initial time to the state of the system is factorized like 
Po = p c ® Pn, where p c is the cavity mode density operator and p n is the 
thermal state of the bath with 

ti{ Pn b\u)b(u')} = n(uj)5{uj - J) = * 8{u - (J), (6.20) 

ex P(^ir) " 1 

now, since the cavity mode has a spectrum centered in uj = u c , within the 
rotating wave approximation we have 



(M*)) = °. ( 6 - 21 ) 

(a\ n {t)a m {t')) ~ n(u c )5(t-t'), (6.22) 
(a in {t)a\ n {t')) ~ [n{u c ) + l]8(t-t') : (6.23) 

moreover, if we deal with optical frequencies, then n(u c ) « 1 and the mean 
excitation number can be usually neglected. 

We observe that the Langevin equation (|6.17p is very general, since it is 
valid not only if the input state is a thermal filed, but e.g. also if it is a 
coherent or squeezed state. Obviously, in the these cases, Eq. (]6. 21116.2216.230 
will be different. 



6.2 Quantum Brownian motion 

The quantum Brownian motion of a massive particle in a potential is more 
intricate than the optical counterpart described in the previous section. The 



Chapter 6 - Quantum Optomechanics 



57 



problem is that the mechanical coupling is different form the optical one and 
moreover we cannot apply the rotating wave approximation; this implies dif- 
ferent properties of the noise operator, which has no more the physical in- 
terpretation of an optical input field but it represents a driving stochastic 
force. 

Let us start with the Hamiltonian of a particle (lower case operators) in a 
bath of mechanical oscillators (upper case operators) , 



H 



P 

2m 



+ V(q) + 



co p/2 t, 



(6.24) 



where the interaction is assumed to depend only on the relative distance 
between the particle and the oscillators: Q'^ — q. 

The following calculations are very similar but not exactly equal to that 
given by Giovannetti and Vitali [15] , because here we start from an equivalent 
but different Hamiltonian (|6.24p . proposed by Gardiner and Zoller [8]. 

Calculations can be simplified using dimensionless quadratures for the bath 
operators, 



so that the Hamiltonian becomes, 
~ 2 1 




1 



duo 



o 



P' 



where the only non- vanishing commutation relations are 



[q,p] 
Qua Pl 



ih , 

i8(u — to') 



(6.25) 
(6.26) 

(6.27) 



(6.28) 
(6.29) 



and = is the coupling constant. 

The Heisenberg equations for the bath operators are 



LOP, 



-ujQ w + K UJ \j-q 



(6.30) 
(6.31) 



and if we define the annihilation operators b w = {Q u + iP^/V^, we have 



-iubu + iK UJ \l —q, 



(6.32) 
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which integrated gives 

MO = e~ M '~* o) Mio)+«u 



(6.33) 



to 



where to is an arbitrary initial time. 

Now we write the Heisenberg equation for the particle variables 



a P 

q = — , 

m 



(6.34) 



P 



- h W,P] 



'fojj 



dujklq+ I dwku\/ —(bu + bl) , (6.35) 



if we substitute (16.330 in (|6.35|) we get 



p(t) = -[v,p(t)]+a(t)+ I 



q(t) +uj l d£ sm[uj{t - t')\q(t') 



to 



where we have defined the Brownian force operator 

' huj r 



dujk UJ \l — 



e -^(*-*o)S a ,(t )+e^(*-*°)St(t ) 



(6.36) 
(6.37) 



The time integral in (|6.36p can be done by parts, giving 



Pit) = -[ V ,p(t)}H(ty 



dujk r . 



cos[a;(t — to)]q(to) + / dt' cos[uj(t — t')] 



to 



m 
(6.38) 

If we suppose a frequency independent coupling = 27 /ir, using that J °° duj cos(. 
irS(t), we get 

kt) = I [V, P(t)} - -P(t) + i(t) - 2 7 5(t - t )q(t ) . (6.39) 

The last term in (|6.39p can be neglected for t > to, therefore we arrive to a 
quantum Langevin equation which is very similar to the classical one 



P(t) = Uv,p(t)]-^p(t)+i(t). 
n m 



(6.40) 



The difference between the QLE (|6.40p and the classical countrepart is on the 
commutation rules and on the second moments of the stochastic force 
(I6.37|) . In fact it can be shown [36] that 



[£(t),i(t')) = 2ih 1 j t 6(t-t') 
and, if the particle is in a thermal bath, then 



«(<)) 



2hj 



7T 



([£(*), f(i')]+> = — - / (ko to cos[u;(t - i')] coth( 



2K B T 



(6.41) 



(6.42) 
), (6.43) 
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It is notable form f (|6.43p . that a truly quantum Brownian motion is not a 
Markov process, however in almost every situation, we can apply the high 
temperature limit 



6.3 Radiation pressure 

J. C. Maxwell, in his Treatise on Electricity and Magnetism of 1873, first 
predicted that light possesses a momentum and therefore it exerts a pressure 
when it hits upon a surface, 

"... in a medium in which waves are propagated, there is a pressure 
in the direction normal to the waves, and numerically equal to the 
energy in unit of volume." [2] 

This pressure is very small, for example the Sun light pressure on Earth is 
of the order of some micropascals. However it becomes extremely important 
in the dynamics of micro/nano optomechanical systems, where the typical 
interactions are very weak and the pressure of a laser field on a mechanical 
device, e.g. on a mirror, is no more negligible. 

We are going to derive the quantum mechanical Hamiltonian of a cavity 
mode coupled with a movable mirror (usually one of the cavity extremes). 



6.3.1 Semiclassical theory of radiation pressure 

We model the system as n photons hitting a movable mirror of an optical 
cavity. 




(6.44) 



and we recover a Markov process 




(6.45) 
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p P h i 


AP m = 2P ph 
> 




i 


*■ 



Oq 



Figura 6.2: Elastic collisions of photons upon the movable mirror Mi- M\ is 
fixed. 



When a photon with momentum p p h = E p h/c makes an elastic collision 
with the mirror (see Fig. I6.2D . the conservation of total momentum imposes a 
change in the mirror momentum of 



Ap m = 2p ph 



2E, 



ph 



(6.46) 



The time delay between different collisions of a photon with the mirror is given 
by the round trip time 



At 



2L 



(6.47) 



The total force exerted by all the n photons on the mirror is given by the 
Newton law 

Aft 



F = n- 



n- 



At L 
corresponding to the classical Hamiltonian 



H — Hught + H v 



+ Hi n t, 



with 



H, 



E, 



int 



-71- 



ph 



L 



(6.48) 

(6.49) 
(6.50) 



where q is the displacement of the mirror form its equilibrium position. 

The Hamiltonian quantization can be done associating to E p h the photon 
energy E p h = huj and to the variables n and q the corresponding quantum 
mechanical operators, 



n 



a) a, 



(6.51) 
(6.52) 
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so that the quantum Hamiltonian becomes 

" — Hugfit + H m i rror + Hi n t, 

with 



int 



fvjj 



a^a q. 



(6.53) 
(6.54) 



6.3.2 Eigenmodes theory of radiation pressure 



Another derivation of the radiation pressure potential (|6.54p can be obtained 
starting from the eigenmodes equation connecting the supported frequencies 
uj n and the length L of the cavity. 




Oq 



Figura 6.3: Variations of the cavity length cause variations of the cavity 
eigenfrequencies . 



In the simple configuration of Fig. (16. 3D , if we suppose the mirrors M\ and 
M% fixed at x = —L and x = respectively, then the standing waves condition 
requires that the eigenmodes are such that nA/2 = L with n £ N. This means 
that the allowed wave vectors are 



2tt n-K 

kn = ~ = T 



(6.55) 



and using the plane wave relation u> n = k n c, we find the eigenmodes equation 

W7TC 



OJr. 



L 



(6.56) 



Let us suppose that we pump a particular mode ui n and that the mirror of the 
cavity are fixed, then the Hamiltonian will be 



(6.57) 



62 



Chapter 6 - Quantum Optomechanics 



Now we relax the condition of a fixed cavity length and we let the right mirror 
Mi free to move around its equilibrium position x = (see Fig. I6.3[) . if we 
call with q « L the displacement of the mirror, then the cavity eigenmodes 
will depend on q through 



mrc mrc 

r^j : 



1 



L 



w fl (0) 



1 



L + q L 

The Hamiltonian of the optomechanical system will be 

H = fmj{q)a) a + H m , 



L 



(6.58) 



(6.59) 



where H m is the Hamiltonian of the non interacting mirror, which can be 
treated for example like a free mass or like a harmonic oscillator, depending 
on the particular system. 

Substituting the linear expansion of the frequency (|6.58p in (|6.59l) , we can 
separate a free field energy and an interaction potential which is exactly equal 
to that given in (|6.54p . 



Cavity with a membrane inside 



— a'aq. 

Li 



(6.60) 



Differently from the semiclassical derivation which is good for perfectly reflec- 
tive mirrors, the expansion of the eigen-frequencies as functions of q can be 
easily generalized to more complex systems, like partially reflective mirrors, 
or cavities with a thick membrane inside [381 [56] . 




q 



Figura 6.4: A thick membrane M3 is placed inside an optical cavity. Radiation 
pressure and vibrations cause small displacements q from the rest position q^. 
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Let us consider a cavity of length 2L and a vibrating membrane positioned 
at x = qo (see Fig. 16. 4p . 

It can be shown that the eigenmodes of the system satisfy the following 
equations 

cot k(L + q ) + cot k(L - q ) = 2y/R/T, (6.61) 

where k is the wave vector, while R and T are the membrane reflectivity and 
transmissivity respectively. 

We first study the special limit qo = T = 0. In this case we have two equal 
cavities and the supported frequencies are the same of (|6,56p . For each integer 
n we have two modes of frequency u n , one on the left and one on the right of 
the membrane. 

Now we relax the condition of perfect reflectivity by considering a mem- 
brane with T 7^ and we assume also a de-centered position of the membrane 
go 7^ 0. The symmetry has been broken, and for each integer n the previous 
eigenmode uj n splits in 



uj n (q ) = + — [sin 1 (V r Rcos2Ai n go) - sin 1 Vr] , (6.62) 
w n(<?o) = oj n + 777 [7T ~ sin~ 1 (v A Rcos2£; n <2o) - sin" 1 Vr]. (6.63) 

These expressions are just a different rewriting of to the implicit relation (|6.61|) . 
Now we consider small displacements q « qo of the membrane around its 
equilibrium position go- Equations (|6.62H6.63]) can be linearized, giving 

^n(q) = u n - 5' - fjr-q, (6-64) 

u+(q) = U ; n + 5 + +f l ^q, (6.65) 



where 



<T = 777 [sin" 1 v A R-sin" 1 (v A Rcos2/c n g )], (6.66) 

5 + = -^-[7r- S m- 1 (VRcos2k n qo) -sin" 1 Vr], (6.67) 
sin2A;„gn 

/ = = (6.68) 

y R 1 — cos 2 2k n qo 

If the equilibrium position go is such that the sine function in (|6.68|) is 

different from zero, then the position dependace of the frequency is very similar 

to (|638jh 

The Hamiltonian of the system, corresponding to the splitted modes (16.641 - 
EMD, is 

H = h(uj n -5-)^a + h(uJn + S + )b Ji b + H m (6.69) 
-nf^a^aq + hf^tfbq, (6.70) 
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where, in the first line there are the free field Hamiltonians of the two shifted 
optical modes (with lowering operators a and b) and of the moving membrane, 
while in the second line there are the radiation pressure contributions of the 
two modes. We observe that, in this configuration, the radiation pressure 
terms have opposite signs corresponding to opposite force directions. 



Capitolo 7 



One driven mode and a 
movable mirror lssi 



7.1 Introduction 

Experimental demonstration of genuine quantum states of macroscopic me- 
chanical resonators with a mass in the nanogram-milligram range represents 
an important step not only for the high-sensitive detection of displacements 
and forces, but also for the foundations of physics. It would represent, in 
fact, a remarkable signature of the quantum behavior of a macroscopic object, 
allowing to shed further light onto the quantum-classical boundary. Signifi- 
cant experimental [371 EH] and theoretical [3D] efforts are currently devoted to 
cooling such microresonators to their quantum ground state. 

However, the generation of other examples of quantum states of a micro- 
mechanical resonator has been also considered recently. The most relevant 
examples are given by squeezed and entangled states. Squeezed states of 
nano-mechanical resonators are potentially useful for surpassing the standard 
quantum limit for position and force detection. 

The conditions under which entanglement between macroscopic objects 
can arise is also currently investigated. After the proposal of Ref. [18], in 
which two mirrors of a ring cavity are entangled by the radiation pressure 
of the cavity mode, other optomechanical systems have been proposed for 
entangling optical and/or mechanical modes by means of the radiation pressure 
interaction. Refs. [33] considered the simplest scheme capable of generating 
stationary optomechanical entanglement, i.e., a single Fabry-Perot cavity with 
a movable mirror. 

Here we shall reconsider the Fabry-Perot model of Ref. [43j, which is re- 
markable for its simplicity and robustness against temperature of the result- 
ing entanglement, and extend its study in various directions. In particular 
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we shall develop a general theory showing how the entanglement between the 
mechanical resonator and optical output modes can be properly defined and 
calculated. This is important since any quantum communication application 
involves traveling output modes rather than intracavity ones, furthermore, by 
considering the output field, one can adopt a multiplexing approach that is, 
using appropriate filters, one can always select many different traveling output 
modes originating from a single intracavity mode (see Fig. 1). 



laser 




Figura 7.1: Scheme of the cavity, which is driven by a laser and has a vibrating 
mirror. With appropriate filters one can select N independent modes from the 
cavity output field. 

We shall see that the relevant dynamics induced by radiation pressure in- 
teraction is carried by the two output modes corresponding to the first Stokes 
and anti-Stokes sidebands of the driving laser. In particular, the optomechan- 
ical entanglement with the intracavity mode is optimally transferred to the 
output Stokes sideband mode, which is however robustly entangled also with 
the anti-Stokes output mode. We shall see that the present Fabry-Perot cavity 
system is preferable with respect to the free space model of Refs. [33], because 
entanglement is achievable in a much more accessible experimental parameter 
region. 

7.2 System dynamics 

We consider a driven optical cavity coupled by radiation pressure to a mi- 
cromechanical oscillator. The typical experimental configuration is a Fabry- 
Perot cavity with one mirror much lighter than the other (see e.g. [38]), but 
our treatment applies to other configurations such as the silica toroidal micro- 
cavity of Refs. [37]. Radiation pressure couples each cavity mode with many 
vibrational normal modes of the movable mirror. However, by choosing the de- 
tection bandwidth so that only an isolated mechanical resonance significantly 
contributes to the detected signal, one can restrict to a single mechanical oscil- 
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lator, since inter-mode coupling due to mechanical nonlinearities are typically 
negligible. The Hamiltonian of the system reads 



H = huj c a)a + -hLo m (p 2 + q 2 ) — hG^a^aq (7.1) 
+ihE{a) e- iwot - ae iuJ0t ). (7.2) 

The first term describes the energy of the cavity mode, with lowering operator 
a ([a, a'] = 1), cavity frequency uj c and decay rate n. The second term gives the 
energy of the mechanical mode, modeled as harmonic oscillator at frequency 
uj m and described by dimensionless position and momentum operators q and 
p ([q,p] = i)- The third term is the radiation-pressure coupling of rate Go = 
(Lu c /L)^/h/muj mi where m is the effective mass of the mechanical mode [37] , 
and L is an effective length that depends upon the cavity geometry: it coincides 
with the cavity length in the Fabry-Perot case, and with the toroid radius 
in the case of Refs. [37] . The last term describes the input driving laser 
with frequency Uq, where E is related to the input laser power P by \E\ = 
^Pk/^o. One can adopt the single cavity mode description of Eq. (17. 2D 
as long as one drives only one cavity mode and the mechanical frequency 
uj m is much smaller than the cavity free spectral range FSR ~ c/L. In this 
case, scattering of photons from the driven mode into other cavity modes is 
negligible [4*8] . 

The dynamic is also determined by the fluctuation-dissipation processes 
affecting both the optical and the mechanical mode. As we have seen in 
the previous sections, these fluctuations can be taken into account using the 
Langevin equations approach (I6.18)6.40p . In this case, we can write a set of 
nonlinear QLE in the interaction picture with respect to ttwoala, 

q = uj m p, (7.3a) 
p = -io m q - 7 m p + Goal a + £, (7.3b) 
a = -(K + iA )a + iG aq + E + V2na in , (7.3c) 

where Ao = iv c — ujq. The mechanical mode is affected by a viscous force with 
damping rate 7 m and by a Brownian stochastic force with zero mean value £, 
that obeys the correlation function (|6.37|) . which rewritten for dimensionless 
mirror operators becomes 



coth (^) +1 



(7.4) 



w m J 2ir 

where ks is the Boltzmann constant and T is the temperature of the reservoir 
of the micromechanical oscillator. The cavity mode amplitude instead decays 
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at the rate k and is affected by the vacuum radiation input noise a in (t), whose 
correlation functions are given by (see l6.2lTl6.23p 

(a in (t)a in *(t')) = [N(lo c ) + 1] 5(t - 0, (7.5) 

and 

(a in ^(t)a in (t , ))=N(u c )6(t-t'), (7.6) 

where N(lo c ) = (exp{Huj c / ksT} — 1) _1 is the equilibrium mean thermal pho- 
ton number. At optical frequencies htOc/ksT 3> 1 and therefore N(uj c ) ~ 0, 
so that only the correlation function of Eq. (|7.5|) is relevant. We shall neglect 
here technical noise sources, such as the amplitude and phase fluctuations of 
the driving laser. They can hinder the achievement of genuine quantum ef- 
fects (see e.g. [17]), but they could be easily accounted for by introducing 
fluctuations of the modulus and of the phase of the driving parameter E of 
Eq. CLl. 

7.2.1 Linearization and stability analysis 

As shown in [43] . significant optomechanical entanglement is achieved when 
radiation pressure coupling is strong, which is realized when the intracavity 
field is very intense, i.e., for high- finesse cavities and enough driving pow- 
er. In this limit (and if the system is stable) the system is characterized by 
a semiclassical steady state with the cavity mode in a coherent state with 
amplitude a s = E/(k + iA), and the micromechanical mirror displaced by 
q s = Go\a s \ 2 /u> m (see Refs. [US 09] for details). The expression giving the in- 
tracavity amplitude a s is actually an implicit nonlinear equation for a s because 

A = A - (7.7) 

is the effective cavity detuning including the effect of the stationary radiation 
pressure. As shown in Refs. [33], when \a s \ 3> 1 the quantum dynamics of 
the fluctuations around the steady state is well described by linearizing the 
nonlinear QLE of Eqs. (17.3|) . Defining the cavity field fluctuation quadratures 
SX = (5a + Sa^)/ y/2 and SY = (5a — 5a^)/iy/2, and the corresponding Hermi- 
tian input noise operators X m = (a in + a in '^) / y2 and Y m = (a m — a %n '^) / 
the linearized QLE can be written in the following compact matrix form [43J 

u(t) = Au(t) + n(t), (7.8) 

where u T (t) = (5q(t), 5p(t), 5X(t), 5Y (t)) T ( T denotes the transposition) is the 
vector of CV fluctuation operators , n T (t) = (0,£(i), y/2KX in (t), y/2KY in (t)) T 
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the corresponding vector of noises and A the matrix 



A 





-U) n 


G 



U). 





7m G 

—k A 

-A -k 



\ 



(7.9) 



where 



G = G a s V2 




Pk 



mw m wo (k 2 + A 2 ) ' 



(7.10) 



is the effective optomechanical coupling (we have chosen the phase reference 
so that a s is real and positive). When a s 3> 1, one has G 3> Go, and therefore 
the generation of significant optomechanical entanglement is facilitated in this 
linearized regime. 

The formal solution of Eq. jj\8]) is u{t) = M(t)u(0) + J Q * dsM(s)n{t - s), 
where M(f) = ex.p{At}. The system is stable and reaches its steady state 
for t — > oo when all the eigenvalues of A have negative real parts so that 
M(oo) = 0. The stability conditions can be derived by applying the Routh- 
Hurwitz criterion [7], yielding the following two nontrivial conditions on the 
system parameters, 

si = 2j m K I k 2 + (uj m - A) 2 k 2 + (uj m + A)' 
+7m [(7m + 2k) {k 2 + A 2 ) + 2/so;^] } 
+Auj m G 2 ( 7m + 2k) 2 > 0, 



s 2 = to m (k 2 + A 2 ) - G 2 A > 0. 



(7.11a) 
(7.11b) 



which will be considered to be satisfied from now on. Notice that when A > 
(laser red-detuned with respect to the cavity) the first condition is always 
satisfied and only s 2 is relevant, while when A < (blue-detuned laser), the 
second condition is always satisfied and only s± matters. 



7.2.2 Correlation matrix of the quantum fluctuations 

The steady state of the bipartite quantum system formed by the vibrational 
mode of interest and the fluctuations of the intracavity mode can be fully char- 
acterized. In fact, the quantum noises £ and a in are zero-mean quantum Gaus- 
sian noises and the dynamics is linearized, and as a consequence, the steady 
state of the system is a zero-mean bipartite Gaussian state, fully characterized 
by its 4 x 4 correlation matrix (CM) VL- = ((uj(oo)%(oo) + Uj(oo)ui(oo))) /2. 
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Starting from Eq. (|7.8p . this steady state CM can be determined in two 
equivalent ways. Using the Fourier transforms Ui(u>) of Ui(t), one has 



Vijit) 



duiduj 
< 

47T 



(u^u^J) + ^(u/RH) . (7.12) 



Then, by Fourier transforming Eq. (17. 8p and the correlation functions of the 
noises, Eqs. (I7,4p and (17. 5|) . one gets 



(ui(w)tij(a/) + Uj(u')ui(u)) 



M{uj)D{u)M{J) t 8(u + u/) 
J ij 



where we have defined the 4x4 matrices 

M[u) = (icu + Ay 1 



(7.13) 
(7.14) 

(7.15) 



and 



D(cj) 



(0 








\ 





7m " coth 

ujm \2k B T J 














K 





\o 








« / 



(7.16) 



The 5(u;+u;') factor is a consequence of the stationarity of the noises, which im- 
plies the stationarity of the CM V: in fact, inserting Eq. (|7.13j) into Eq. (|7.12l) . 
one gets that V is time-independent and can be written as 



V 



dwM{u)D{u)M{io)^ . 



(7.17) 



It is however reasonable to simplify this exact expression for the steady state 
CM, by appropriately approximating the thermal noise contribution I^^) 
in Eq. (|7. 16j) . In fact ksT/h ~ 10 11 s _1 even at cryogenic temperatures and 
it is therefore much larger than all the other typical frequency scales, which 
are at most of the order of 10 9 Hz. The integrand in Eq. (|7.17p goes rapidly 
to zero at uj ~ 10 11 Hz, and therefore one can safely neglect the frequency 
dependence of L>22(w) by approximating it with its zero- frequency value 



7mW 



coth 



\2k B T) 



In 



2k B T 



7 m (2n + 1) 



(7.18) 



where n = (ex.p{Hio m /kBT} — 1) 1 is the mean thermal excitation number of 
the resonator. 
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It is easy to verify that assuming a frequency-independent diffusion ma- 
trix D is equivalent to make the following Markovian approximation on the 
quantum Brownian noise £(t), 

(mm + £(*')£(*)} /2 * lm(2n + l)8(t - t'). (7.19) 

Within this Markovian approximation, the above frequency domain treatment 
is equivalent to the time domain derivation considered in [43] which, starting 
from the formal solution of Eq. (17. 8h , arrives at 

rco /■oo 

V i3 (oo) = y2 ds ds'Mik^M^Dkiis-s'), (7.20) 

where Dki(s-s') = ((wfc(s)raj(s') + ni(s')nk(s)}) /2 is the matrix of the station- 
ary noise correlation functions. The Markovian approximation of the thermal 
noise on the mechanical resonator yields D^i{s — s') = Dki5{s — s'), with 
D = Diag[0, 7 m (2n + 1), k, k), so that Eq. (|7.20|) becomes 

/"OO 

V= dsM(s)DM(s) T , (7.21) 
Jo 

which is equivalent to Eq. (|7.17p whenever D does not depend upon cu. When 
the stability conditions are satisfied (M(oo) = 0), Eq. (|7.2ip is equivalent to 
the following Lyapunov equation for the steady-state CM, 

AV + VA T = -D, (7.22) 

which is a linear equation for V and can be straightforwardly solved, but the 
general exact expression is too cumbersome and will not be reported here. 



7.3 Optomechanical entanglement with the intra- 
cavity mode 

In order to establish the conditions under which the optical mode and the 
mirror vibrational mode are entangled we consider the logarithmic negativity 
Ej\f, which we have defined in (|3.26|) as 

Em = max[0,-ln277~], (7.23) 

where rf = 2~ 1 / 2 E(V) - [S(F) 2 - 4detF] 1/2 , with E(V) = detV m + 
det V c — 2 det V mc , and we have used the 2x2 block form of the CM 

vJfi M. (7.24) 
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7.3.1 Correspondence with the down-conversion process 



As already shown in |37} HO] many features of the radiation pressure interaction 
in the cavity can be understood by considering that the driving laser light is 
scattered by the vibrating cavity boundary mostly at the first Stokes (uQ — u> m ) 
and anti-Stokes (ujq + <*> m ) sidebands. Therefore we expect that the optome- 
chanical interaction and eventually entanglement will be enhanced when the 
cavity is resonant with one of the two sidebands, i.e., when A = =Lu; m . 

It is useful to introduce the mechanical annihilation operator 5b = (5q + 
iSp) I y/2, obeying the following QLE 



5b = -iuj m 5b -^(Sb- 5b^ + i- (Sat + 5a) + (7.25) 

Moving to another interaction picture by introducing the slowly-moving tilded 
operators 5b(t) = 5b(t)e~ luJmt and 5a(t) = 5a(t)e~ lAt , we obtain from the 
linearized version of Eq. (|7.3b ) and Eq. (17.250 the following QLEs 



5b = -J^^Sb-SPe^j+^b 



in 



7m ^xi c/.t />U;„ l'*> 

+£ (SaU i ( A +*J t + Sae^™"^ (7.26) 
5h = -K5h + £ (»t e i(A+«m)* + S l e i(A-u m )t^ 

+V2^a in . (7.27) 

Note that we have introduced two noise operators: i) a m (t) = a m (i)e* A * , 
possessing the same correlation function as a m (t); ii) b tn {t) = £(t)e luJmt / 'y/2 
which, in the limit of large u m , acquires the correlation functions [8] 

(b in '\t)b in (t')) = n5(t-t'), (7.28) 
(b in (t)b ir ^(t')} = [n + l]5(t-t'). (7.29) 

Eqs. (fT26]l - (f7T271) are still equivalent to the linearized QLEs of Eq. (fTSjl . but 
now we particularize them by choosing A = ±w m . If the cavity is resonant 
with the Stokes sideband of the driving laser, A = — u m , one gets 

Sb = -— 56 + —5$e 2iuj,nt + i—5a [ 
2 2 2 

C 



+ i-5~ae ZlUmt + ^b in , (7.30) 
5h = - K 5ti + i^5tf + i^5be 2iuJmt + V2^a in , (7.31) 
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while when the cavity is resonant with the anti-Stokes sideband of the driving 
laser, A = u m , one gets 

Si = -lULSb + ^tfe^ + i-Sa 
2 2 2 

+ z^We- 2 ^™* + ^b in , (7.32) 
5h = -K5a + i^5b + i^5tfe- 2iuJmt + V2^a in . (7.33) 

From Eqs. (|7,30p - (|7.3ip we see that, for a blue-detuned driving laser, A = 
— u m , the cavity mode and mechanical resonator are coupled via two kinds 
of interactions: i) a down-conversion process characterized by 5b>5a) + 5a5b, 
which is resonant and ii) a beam-splitter-like process characterized by 5¥5a + 
5a* 5b, which is off resonant. Since the beam splitter interaction is not able to 
entangle modes starting from classical input states, and it is also off-resonant 
in this case, one can invoke the rotating wave approximation (RWA) (which 
is justified in the limit of uj m S> G, k) and simplify the interaction to a down 
conversion process, which is known to generate bipartite entanglement. In the 
red-detuned driving laser case, Eqs. (|7.32p - (|7.33|) show that the two modes are 
strongly coupled by a beam-splitter-like interaction, while the down-conversion 
process is off-resonant. If one chose to make the RWA in this case, one would be 
left with an effective beam splitter interaction which cannot entangle. There- 
fore, in the RWA limit uj m 3> G, k, the best regime for strong optomechanical 
entanglement is when the laser is blue-detuned from the cavity resonance and 
down-conversion is enhanced. However, as it will be seen in the following sec- 
tion, this is hindered by instability and one is rather forced to work in the 
opposite regime of a red-detuned laser where instability takes place only at 
large values of G. 



7.3.2 Entanglement in the blue-detuned regime 



The CM of the Gaussian steady state of the bipartite system, can be obtained 
from Eqs. (|730j) - (|73ll and Eqs. Eqs. (17321 - (17331) in the RWA limit, with the 
techniques of the former section (see also |41| ) 



v = v- 



( 












\ 







V?i 
















vi 







\ 










^33 


/ 



(7.34) 
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where the upper (lower) sign corresponds to the blue- (red-)detuned case, and 



11 2 ( 7m + 2 K )(2 7m At T G 2 )' v ; 

Vi = -+^ [ ^ 1/2±1/ 1, (7.35b) 
33 2 ( 7m + 2 K )(2 7m ^G' 2 )' 1 7 

± 2G 7m K[n + l/2±l/2] 

14 ( 7m + 2K)(2 7mKT G2)' { - ) 

For clarity we have included the red-detuned case in the RWA approximation 
and we see that det V^ c = TO^u) 2 ) i- e -i is non-negative in this latter case, 
which is a sufficient condition for the separability of bipartite states [16 1. Of 
course, this is expected, since it is just the beam-splitter interaction that 
generates this CM. Thus, in the weak optomechanical coupling regime of the 
RWA limit, entanglement is obtained only for a blue-detuned laser, A = —ui m . 
However, the amount of achievable optomechanical entanglement at the steady 
state is seriously limited by the stability condition of Eq. (|7.11a ). which in the 



RWA limit A = — u m 3> K,J m , simplifies to G < y/2Kj m . Since one needs 
small mechanical dissipation rate 7m in order to see quantum effects, this 
means a very low maximum value for G. The logarithmic negativity Ej\f is an 
increasing function of the effective optomechanical coupling G (as expected) 
and therefore the stability condition puts a strong upper bound also on Ej\f. 
In fact, it is possible to prove that the following bound on Ej\f exists 



Em < In 



1 + n 



(7.36) 



showing that Ej\f < In 2 and above all that entanglement is extremely fragile 
with respect to temperature in the RWA limit because, due to the stability 
condition, Ej^ vanishes as soon as n > 1. 



7.3.3 Entanglement in the red-detuned regime 

We conclude that, due to instability, one can find significant optomechani- 
cal entanglement, which is also robust against temperature, only far from the 
RWA regime, in the strong coupling regime in the region with positive A, be- 
cause Eq. (|7.11bj) allows for higher values of G. This is confirmed by Fig. 17.21 
where Ej\f is plotted versus the normalized detuning A/u m and the normal- 
ized input power P/Po, (Pq = 50 mW) at a fixed value of the cavity finesse 
F = F = 1.67 x 10 4 in (a), and versus the normalized finesse F/Fq and nor- 
malized input power P/Pq at a fixed cavity detuning A = to m in (b). We have 
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assumed an experimentally achievable situation, i.e., a mechanical mode with 
uj m /2TT = 10 MHz, Q = 10 5 , mass m = 10 ng, and a cavity of length L = 1 
mm, driven by a laser with wavelength 810 nm, yielding Go = 0.95 kHz and 
a cavity bandwidth k = 0.9uj m when F = Fq. We have assumed a reservoir 
temperature for the mirror T = 0.4 K, corresponding to n ~ 833. Fig. 17.2a 
shows that Ej^ is peaked around A ~ u m , even though the peak shifts to 
larger values of A at larger input powers P. For increasing P at fixed A, 
Ej^ increases, even though at the same time the instability region (where the 
plot suddenly interrupts) widens. In Fig. 17.2b we have fixed the detuning at 
A = ijj m (i.e., the cavity is resonant with the anti-Stokes sideband of the laser) 
and varied both the input power and the cavity finesse. We see again that 
Ej^f is maximum just at the instability threshold and also that, once that the 
finesse has reached a sufficiently large value, F ~ Fq, Ej\f roughly saturates at 
larger values of F. That is, one gets an optimal optomechanical entanglement 
when k ~ <jj m and moving into the well-resolved sideband limit n <C co m does 
not improve the value of Ej\f. The parameter region analyzed is analogous 
to that considered in [33], where it has been shown that this optomechani- 
cal entanglement is extremely robust with respect to the temperature of the 
reservoir of the mirror, since it persists to more than 20 K. 
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Figura 7.2: (Color online) (a) Logarithmic negativity Ej^ versus the normal- 
ized detuning A/uj m and normalized input power P/Pq, (Pq = 50 mW) at 
a fixed value of the cavity finesse F = Fq = 1.67 x 10 4 ; (b) Ej\f versus the 
normalized finesse F/Fq and normalized input power P/Pq at a fixed detuning 
A = uj m . Parameter values are u m /27r = 10 MHz, Q = 10 5 , mass m = 10 
ng, a cavity of length L = 1 mm driven by a laser with wavelength 810 nm, 
yielding Go = 0.95 KHz and a cavity bandwidth k = 0.9w m when F = Fq. We 
have assumed a reservoir temperature for the mirror T = 0.4 K, correspond- 
ing to n ~ 833. The sudden drop to zero of Ej^f corresponds to entering the 
instability region. 
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7.3.4 Relationship between entanglement and cooling 

As discussed in detail in [3D] the same cavity-mechanical resonator system 
can be used for realizing cavity-mediated optical cooling of the mechanical 
resonator via the back-action of the cavity mode. In particular, back-action 
cooling is optimized just in the same regime where A ~ uj m . This fact is 
easily explained by taking into account the scattering of the laser light by the 
oscillating mirror into the Stokes and anti-Stokes sidebands. The generation of 
an anti-Stokes photon takes away a vibrational phonon and is responsible for 
cooling, while the generation of a Stokes photon heats the mirror by producing 
an extra phonon. If the cavity is resonant with the anti-Stokes sideband, 
cooling prevails and one has a positive net laser cooling rate given by the 
difference of the scattering rates. 

It is therefore interesting to discuss the relation between optimal optome- 
chanical entanglement and optimal cooling of the mechanical resonator. This 
can easily performed because the steady state CM V determines also the res- 
onator energy, since the effective stationary excitation number of the resonator 
is given by n e // = (Vii + V22 — 1) /2 (see Ref. |30] for the exact expression of 
these matrix elements giving the steady state position and momentum res- 
onator variances). In Fig. 17.31 we have plotted n e tt under exactly the same pa- 
rameter conditions of Fig. 17.21 We see that ground state cooling is approached 
(n e ff < 1) simultaneously with a significant entanglement. This shows that 
a significant back-action cooling of the resonator by the cavity mode is an 
important condition for achieving an entangled steady state which is robust 
against the effects of the resonator thermal bath. 

Nonetheless, entanglement and cooling are different phenomena and opti- 
mizing one does not generally optimize also the other. This can be seen by 
comparing Figs. 17.21 and 17. 3t Ej\f is maximized always just at the instability 
threshold, i.e., for the maximum possible optomechanical coupling, while this 
is not true for n e ff, which is instead minimized quite far from the instability 
threshold. For a more clear understanding we make use of some of the re- 
sults obtained for ground state cooling in Refs. [ID]. In the perturbative limit 
where G <C uj m ,K, one can define scattering rates into the Stokes (A + ) and 
anti-Stokes (A-) sidebands as 



A± 



G 2 k/2 



(7.37) 




so that the net laser cooling rate is given by 



r = A 



A+ > 0. 



(7.38) 
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The final occupancy of the mirror mode is consequently given by [3D] 



im i L I m T L 



where the first term in the right hand side of the above equation is the 
minimized thermal noise, that can be made vanishingly small provided that 
7m "C I\ while the second term shows residual heating produced by Stokes 
scattering off the vibrational ground state. When V ^> ^mfi, the lower bound 
for n e ff is practically set by the ratio A + /T. However, as soon as G is in- 
creased for improving the entanglement generation, scattering into higher or- 
der sidebands takes place, with rates proportional to higher powers of G. As 
a consequence, even though the effective thermal noise is still close to zero, 
residual scattering off the ground state takes place at a rate that can be much 
higher than A+. This can be seen more clearly in the exact expression of 
(5q 2 ) = Vn given in |40j, which is shown to diverge at the threshold given by 
Eq. $njB) . 
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Figura 7.3: (Color online) (a) Effective stationary excitation number of the 
resonator n e ff versus the normalized detuning A/uj m and normalized input 
power P/Pq, (Po = 50 mW) at a fixed value of the cavity finesse F = Fq = 
1.67 x 10 4 ; (b) n e fj versus the normalized finesse F/Fq and normalized input 
power P/Pq at a fixed detuning A = w m . Parameter values are the same 
as in Fig. 17.21 Again, the sudden drop to zero corresponds to entering the 
instability region. 



7.4 Optomechanical entanglement with cavity out- 
put modes 

The above analysis of the entanglement between the mechanical mode of in- 
terest and the intracavity mode provides a detailed description of the internal 
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dynamics of the system, but it is not of direct use for practical applications. 
In fact, one typically does not have direct access to the intracavity field, but 
one detects and manipulates only the cavity output field. For example, for any 
quantum communication application, it is much more important to analyze 
the entanglement of the mechanical mode with the optical cavity output, i.e., 
how the intracavity entanglement is transferred to the output field. Moreover, 
considering the output field provides further options. In fact, by means of 
spectral filters, one can always select many different traveling output modes 
originating from a single intracavity mode and this gives the opportunity to 
easily produce and manipulate a multipartite system, eventually possessing 
multipartite entanglement. 

7.4.1 General definition of cavity output modes 

The intracavity field Sa(t) and its output are related by the input-output 
relation given in (|6.19|) . which rewritten with the notation of this chapter 
becomes 

a ou \t) = V2^8a(t) - a in (t), (7.40) 

where the output field possesses the same correlation functions of the optical 
input field a tn (t) and the same commutation relation, i.e., the only nonzero 
commutator is [a out (t), a out (t')^ =5(t — t'). From the continuous output field 
a out (t) one can extract many independent optical modes, by selecting different 
time intervals or equivalently, different frequency intervals (see e.g. |50|). One 
can define a generic set of N output modes by means of the corresponding 
annihilation operators 

a° k u \t)= f dsg k (t-s)a out (s), k = l,...N, (7.41) 

J — oo 

where gk(s) is the causal filter function defining the k-th output mode. These 
annihilation operators describe N independent optical modes when a?"*(i), a k mt (t) j; 
Sjk, which is verified when 







dsgj{s)*g k (s) = S jk , (7.42) 



i.e., the N filter functions gk(t) form an orthonormal set of square-integrable 
functions in [0, oo). The situation can be equivalently described in the fre- 
quency domain: taking the Fourier transform of Eq. (|7.4ip . one has 

/oo JJ. 
-oo V^TT 



Chapter 7 - One driven mode and a movable mirror 



81 



where gk(w) is the Fourier transform of the filter function. An explicit example 
of an orthonormal set of filter functions is given by 

9k® = W-W-rK -*** (7.44) 

(9 denotes the Heavyside step function) provided that f2fc and r satisfy the 
condition 

Qj — £lk = — p, integer p. (7-45) 

These functions describe a set of independent optical modes, each centered 
around the frequency Ofe and with time duration r, i.e., frequency bandwidth 
~ 1/t, since 

JL J(u-n k )r/2 s ' m [(w - fifc)r/2] 



When the central frequencies differ by an integer multiple of 2tt/t, the cor- 
responding modes are independent due to the destructive interference of the 
oscillating parts of the spectrum. 

7.4.2 Stationary correlation matrix of output modes 

The entanglement between the output modes defined above and the mechan- 
ical mode is fully determined by the corresponding (2N + 2) x (2N + 2) CM, 
which is defined by 

V**(t) = \ {u™\t)uf\t) + uf \t)uT\t)) , (7.47) 

where 

u out (t) (7.48) 

= (j 3 (t),^),irw, Yr\t), x° N u \t),Y^\t)) T 

is the vector formed by the mechanical position and momentum fluctua- 
tions and by the amplitude {X° k ut {t) = [a° k ut {t) + a° k ut [tf] fy/2), and phase 
(Y k out {t) = [a° k ut (t) - a° k ut {tf] /iy/2) quadratures of the N output modes. The 
vector -u olt *(t) properly describes N + 1 independent CV bosonic modes, and in 
particular the mechanical resonator is independent of (i.e., it commutes with) 
the N optical output modes because the latter depend upon the output field 
at former times only (s < t). From the definition of u out (t), of the output 
modes of Eq. (|7.41|) , and the input-output relation of Eq. (|7.4U|) one can write 

v™\t) = f dsT ik (t - s)u% xt {s) 



dsT ik {t - s)n e k x \ S ), (7.49) 

oo 
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where 

u ext (t) = (6q(t),Sp(t),X(t), Y(t),..., X(t), Y(t)) T (7.50) 

is the 2N + 2-dimensional vector obtained by extending the four-dimensional 
vector u(t) of the preceding section by repeating N times the components 
related to the optical cavity mode, and 

1 



, ext 



n 



(t) 



(0, 0, X in (t), Y in (t), . . . , X in (t), Y in {t)) J 



(7.51) 



is the analogous extension of the noise vector n{t) of the former section without 
however the noise acting on the mechanical mode. In Eq. (|7.49p we have 
also introduced the (2N + 2) x (2N + 2) block-matrix consisting of N + 1 
two-dimensional blocks 



T(t) = 



( S(t) 
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S(t) 
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Mmpi(i) 
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2Klm<72(0 
2nReg 2 (t) 



(7.52) 

Using Fourier transforms, and the correlation function of the noises, one can 
derive the following general expression for the stationary output correlation 
matrix, which is the counterpart of the 4x4 intracavity relation of Eq. (|7.17|) 



V 



out 



dojT{uj) 



M ext {uj) + 



P< 



out 



2k 



D ex \u) 



P, 



out 



2k 



T( W )t, 
(7.53) 

where P ou t = Diag[0, 0, 1,1,.. .] is the projector onto the 2A r -dimensional space 
associated with the output quadratures, and we have introduced the extensions 
corresponding to the matrices M{oj) and D (u) of the former section, 



M ext {uo) = (iuj + A ext ) 



(7.54) 



with 
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and 



D ex \u) 



/o 







V '■ 



1 m 



coth 
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(7.56) 



/ 



A deeper understanding of the general expression for V out of Eq. (|7.53|) is 
obtained by multiplying the terms in the integral: one gets 



V 



out 



dujf{uj)M ext {uj)D ext {uj)M ext {uj) ] f{uj) ] + — 



+ 



- 2 I cIujT(uj) M ext (u;)R ou t + Roui- 



M ex \u)^ f(u)\ (7.57) 



where R ou t = PoutD (uj)/k = D {ui)P out / k and we have used the fact that 



/ 



^f(u;)P out D £Xt (u,)P out f(^ = ?f-. 



(7.58) 



The first integral term in Eq. (17.570 is the contribution coming from the in- 
teraction between the mechanical resonator and the intracavity field. The 
second term gives the noise added by the optical input noise to each output 
mode. The third term gives the contribution of the correlations between the 
intracavity mode and the optical input field, which may cancel the destructive 
effects of the second noise term and eventually, even increase the optomechan- 
ical entanglement with respect to the intracavity case. We shall analyze this 
fact in the following section. 



7.4.3 A single output mode 

Let us first consider the case when we select and detect only one mode at the 
cavity output. Just to fix the ideas, we choose the mode specified by the filter 
function of Eqs. (I7.44p and (|7.46p . with central frequency Q and bandwidth 
r -1 . Straightforward choices for this output mode are a mode centered either 
at the cavity frequency, O, = oj c — Uq, or at the driving laser frequency, O = 
(we are in the rotating frame and therefore all frequencies are referred to the 
laser frequency loq), and with a bandwidth of the order of the cavity band- 
width t -1 ~ k. However, as discussed above, the motion of the mechanical 
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resonator generates Stokes and anti-Stokes motional sidebands, consequently 
modifying the cavity output spectrum. Therefore it may be nontrivial to de- 
termine which is the optimal frequency bandwidth of the output field which 
carries most of the optomechanical entanglement generated within the cavity. 
The cavity output spectrum associated with the photon number fluctuations 
S(oj) = (Sa(uy Sa(u)) is shown in Fig. 17.41 where we have considered a param- 
eter regime close to that considered for the intracavity case, i.e., an oscillator 
with uj rn /2Tr = 10 MHz, Q = 10 5 , mass m = 50 ng, a cavity of length L = 1 mm 
with finesse F = 2 x 10 4 , detuning A = u m , driven by a laser with input power 
P = 30 mW and wavelength 810 nm, yielding G = 0.43 kHz, G = 0.41cj m , 
and a cavity bandwidth k = 0.75u; m . We have again assumed a reservoir tem- 
perature for the mirror T = 0.4 K, corresponding to n — 833. This regime 
is not far but does not corresponds to the best intracavity optomechanical 
entanglement regime discussed in Sec. 17.31 In fact, optomechanical entangle- 
ment monotonically increases with the coupling G and is maximum just at 
the bistability threshold, which however is not a convenient operating point. 
We have chosen instead a smaller input power and a larger mass, implying a 
smaller value of G and an operating point not too close to threshold. 
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Figura 7.4: (Color online) Cavity output spectrum in the case of an oscillator 
with uj m /2 , K = 10 MHz, Q = 10 5 , mass m = 50 ng, a cavity of length L = 1 
mm with finesse F = 2 x 10 4 , detuning A = u> m , driven by a laser with 
input power P = 30 mW and wavelength 810 nm, yielding Go = 0.43 kHz, 
G = 0.41u; m , and a cavity bandwidth k = 0.75o> TO . We have again assumed a 
reservoir temperature for the mirror T = 0.4 K, corresponding to n ~ 833. In 
this regime photons are scattered only at the two first motional sidebands, at 
wo ± w m . 
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In order to determine the output optical mode which is better entangled 
with the mechanical resonator, we study the logarithmic negativity E_\r asso- 
ciated with the output CM V out of Eq. ([737]> (for N = 1) as a function of 
the central frequency of the mode Q and its bandwidth r _1 , at the consid- 
ered operating point. The results are shown in Fig. 17. b\ where Ej\f is plotted 
versus Q/ui m at five different values of e = Tu> m . If e < 1, i.e., the band- 
width of the detected mode is larger than u m , the detector does not resolve 
the motional sidebands, and Ej^ has a value (roughly equal to that of the in- 
tracavity case) which does not essentially depend upon the central frequency. 
For smaller bandwidths (larger e), the sidebands are resolved by the detec- 
tion and the role of the central frequency becomes important. In particular 
Ej\f becomes highly peaked around the Stokes sideband Q = — u m , showing 
that the optomechanical entanglement generated within the cavity is mostly 
carried by this lower frequency sideband. What is relevant is that the op- 
tomechanical entanglement of the output mode is significantly larger than its 
intracavity counterpart and achieves its maximum value at the optimal value 
e ~ 10, i.e., a detection bandwidth r _1 ~ cj m /10. This means that in practice, 
by appropriately filtering the output light, one realizes an effective entangle- 
ment distillation because the selected output mode is more entangled with the 
mechanical resonator than the intracavity field. 
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Figura 7.5: (Color online) Logarithmic negativity Ej^ of the CV bipartite 
system formed by mechanical mode and a single cavity output mode versus 
the central frequency of the detected output mode £l/uj m at five different 
values of its inverse bandwidth e = uj m T. The other parameters are the same 
as in Fig. 17.41 When the bandwidth is not too large, the mechanical mode is 
significantly entangled only with the first Stokes sideband at loq — uo m . 

The fact that the output mode which is most entangled with the mechani- 
cal resonator is the one centered around the Stokes sideband is also consistent 
with the physics of two previous models analyzed in Refs. [5 1 (, 144] . In [51] an 
atomic ensemble is inserted within the Fabry-Perot cavity studied here, and 
one gets a system showing robust tripartite (atom-mirror-cavity) entanglement 
at the steady state only when the atoms are resonant with the Stokes sideband 
of the laser. In particular, the atomic ensemble and the mechanical resonator 
become entangled under this resonance condition, and this is possible only if 
entanglement is carried by the Stokes sideband because the two parties are 
only indirectly coupled through the cavity mode. In [33] , a free-space optome- 
chanical model is discussed, where the entanglement between a vibrational 
mode of a perfectly reflecting micro-mirror and the two first motional side- 
bands of an intense laser beam shined on the mirror is analyzed. Also in that 
case, the mechanical mode is entangled only with the Stokes mode and it is 
not entangled with the anti-Stokes sideband. 

By looking at the output spectrum of Fig. 17.41 one can also understand 
why the output mode optimally entangled with the mechanical mode has a 
finite bandwidth r _1 ~ u; m /10 (for the chosen operating point). In fact, the 
optimal situation is achieved when the detected output mode overlaps as best 
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as possible with the Stokes peak in the spectrum, and therefore r _1 coincides 
with the width of the Stokes peak. This width is determined by the effective 
damping rate of the mechanical resonator, 7™ ^ = 7m + I\ given by the sum 
of the intrinsic damping rate "y m and the net laser cooling rate T of Eq. (17.381 
It is possible to check that, with the chosen parameter values, the condition 
e = 10 corresponds to r~ 



7; 



It is finally important to analyze the robustness of the present optome- 
chanical entanglement with respect to temperature. As discussed above and 
shown in [43] . the entanglement of the resonator with the intracavity mode is 
very robust. It is important to see if this robustness is kept also by the optome- 
chanical entanglement of the output mode. This is shown by Fig. 17.61 where 
the entanglement Ej\f of the output mode centered at the Stokes sideband 
= — ijj m is plotted versus the temperature of the reservoir at two different 
values of the bandwidth, the optimal one e = 10, and at a larger bandwidth 
e = 2. We see the expected decay of Ej\f for increasing temperature, but 
above all that also this output optomechanical entanglement is robust against 
temperature because it persists even above liquid He temperatures, at least in 
the case of the optimal detection bandwidth e = 10. 
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Figura 7.6: (Color online) Logarithmic negativity Ej^ of the CV bipartite 
system formed by mechanical mode and the cavity output mode centered 
around the Stokes sideband f2 = — Lo m versus temperature for two different 
values of its inverse bandwidth e = Lo m T. The other parameters are the same 
as in Fig. 17.41 
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7.4.4 Two output modes 

Let us now consider the case where we detect at the output two independent, 
well resolved, optical output modes. We use again the step-like filter functions 
of Eqs. (|7.44p and (|7.46j) . assuming the same bandwidth r" 1 for both modes 
and two different central frequencies, f^i and satisfying the orthogonality 
condition of Eq. (|7.45jl Oi — O2 = 2pnT~ 1 for some integer p, in order to have 
two independent optical modes. It is interesting to analyze the stationary 
state of the resulting tripartite CV system formed by the two output modes 
and the mechanical mode, in order to see if and when it is able to show i) 
purely optical bipartite entanglement between the two output modes; ii) fully 
tripartite optomechanical entanglement. 

The generation of two entangled light beams by means of the radiation 
pressure interaction of these fields with a mechanical element has been al- 
ready considered in various configurations. In Ref. [52], and more recently in 
Ref. [42|, two modes of a Fabry- Perot cavity system with a movable mirror, 
each driven by an intense laser, are entangled at the output due to their com- 
mon ponderomotive interaction with the movable mirror. In the single mirror 
free-space model of Ref. [44] , the two first motional sidebands are also robust- 
ly entangled by the radiation pressure interaction as in a two-mode squeezed 
state produced by a non-degenerate parametric amplifier. 

The situation considered here is significantly different from that of Refs. [52} 
I42j . which require many driven cavity modes, each associated with the cor- 
responding output mode. In the present case instead, the different output 
modes originate from the same single driven cavity mode, and therefore it is 
much simpler from an experimental point of view. The present scheme can 
be considered as a sort of "cavity version" of the free-space case of Ref. |44| . 
where the reflecting mirror is driven by a single intense laser. Therefore, as 
in [31], one expects to find a parameter region where the two output modes 
centered around the two motional sidebands of the laser are entangled. This 
expectation is clearly confirmed by Fig. 17.71 where the logarithmic negativity 
Ej\f associated with the bipartite system formed by the output mode centered 
at the Stokes sideband = —oo m ) and a second output mode with the same 
inverse bandwidth (e = w m r = 107r) and a variable central frequency $7, is 
plotted versus Q/u; TO . Ej\f is calculated from the CM of Eq. (|7.57p (for N = 2), 
eliminating the first two rows associated with the mechanical mode, and as- 
suming the same parameters considered in the former subsection for the single 
output mode case. One can clearly see that bipartite entanglement between 



Chapter 7 - One driven mode and a movable mirror 



89 



the two cavity outputs exists only in a narrow frequency interval around the 
anti-Stokes sideband, f2 = u m , where Ej/ achieves its maximum. This shows 
that, as in [44j . the two cavity output modes corresponding to the Stokes and 
anti-Stokes sidebands of the driving laser are significantly entangled by their 
common interaction with the mechanical resonator. The advantage of the 
present cavity scheme with respect to the free-space case of [H] is that the pa- 
rameter regime for reaching radiation-pressure mediated optical entanglement 
is much more promising from an experimental point of view because it requires 
less input power and a not too large mechanical quality factor of the resonator. 
In Fig. 17.81 the dependence of E/j- of the two output modes centered at the 
two sidebands Q = ±tu m upon their inverse bandwidth e is studied. We see 
that, differently from optomechanical entanglement of the former subsection, 
the logarithmic negativity of the two sidebands always increases for decreas- 
ing bandwidth, and it achieves a significant value (~ 1), comparable to that 
achievable with parametric oscillators, for very narrow bandwidths. This fact 
can be understood from the fact that quantum correlations between the two 
sidebands are established by the coherent scattering of the cavity photons by 
the oscillator, and that the quantum coherence between the two scattering 
processes is maximal for output photons with frequencies ujq ± oj m . 
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Figura 7.7: (Color online) Logarithmic negativity Ej^ of the bipartite system 
formed by the output mode centered at the Stokes sideband (Qi = — uj m ) and 
a second output mode with the same inverse bandwidth {e = Lo m T = Wir) and 
a variable central frequency £1, plotted versus £l/uj m . The other parameters 
are the same as in Fig. 17.41 Optical entanglement is present only when the 
second output mode overlaps with the anti-Stokes sideband. 
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Figura 7.8: (Color online) Logarithmic negativity E^f of the bipartite sys- 
tem formed by the two output modes centered at the Stokes and anti-Stokes 
sidebands (0 = ±u> m ) versus the inverse bandwidth e = uj m T. The other 
parameters are the same as in Fig. 17.41 




In Fig. 17.91 we analyze the robustness of the entanglement between the 
Stokes and anti-Stokes sidebands with respect to the temperature of the me- 
chanical resonator, by plotting, for the same parameter regime of Fig. 17.81 
Ej\f versus the temperature T at two different values of the inverse bandwidth 
(s = 107r, IOOtt). We see that this purely optical CV entanglement is ex- 
tremely robust against temperature, especially in the limit of small detection 
bandwidth, showing that the effective coupling provided by radiation pressure 
can be strong enough to render optomechanical devices with high-quality res- 
onator a possible alternative to parametric oscillators for the generation of 
entangled light beams for CV quantum communication. 
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Figura 7.9: (Color online) Logarithmic negativity Ej\f of the two output modes 
centered at the Stokes and anti-Stokes sidebands (f2 = ±uj m ) versus the 
temperature of the resonator reservoir, at two different values of the inverse 
bandwidth e = uj m T. The other parameters are the same as in Fig. 17.41 



Since in Figs. 17.71 and 17.81 we have used the same parameter values for 
the cavity-resonator system used in Fig. 17.51 we have that in this parameter 
regime, the output mode centered around the Stokes sideband mode shows 
bipartite entanglement simultaneously with the mechanical mode and with 
the anti-Stokes sideband mode. This fact suggests that, in this parameter 
region, the CV tripartite system formed by the output Stokes and anti-Stokes 
sidebands and the mechanical resonator mode could be characterized by a fully 
tripartite-entangled stationary state. This is confirmed by Fig. 17.101 where we 
have applied the classification criterion of Ref. [53] , providing a necessary and 
sufficient criterion for the determination of the entanglement class in the case 
of tripartite CV Gaussian states, which is directly computable in terms of the 
eigenvalues of appropriate test matrices [53]. These eigenvalues are plotted 
in Fig. 17.101 versus the inverse bandwidth e at A = u) m in the left plot, and 
versus the cavity detuning A/uj m at the fixed inverse bandwidth e = ir in the 
right plot (the other parameters are again those of Fig. I7.4p . We see that 
all the eigenvalues are negative in a wide interval of detunings and detection 
bandwidth of the output modes, showing, as expected, that we have a fully 
tripartite-entangled steady state. 
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Figura 7.10: (Color online) Analysis of tripartite entanglement. The minimum 
eigenvalues after partial transposition with respect to the Stokes mode (blue 
line), anti-Stokes mode (green line) and the mechanical mode (red line) are 
plotted versus the inverse bandwidth e at A = oj m in the left plot, and versus 
the cavity detuning A/uj m at the fixed inverse bandwidth e = it in the right 
plot. The other parameters are the same as in Fig. 17.41 These eigenvalues 
are all negative in the studied intervals, showing that one has fully tripartite- 
entanglement. 

Therefore, if we consider the system formed by the two cavity output fields 
centered around the two motional sidebands at loq ± uj m and the mechanical 
resonator, we find that the entanglement properties of its steady state are 
identical to those of the analogous tripartite optomechanical free-space system 
of Ref. [44] . In fact, the Stokes output mode shows bipartite entanglement both 
with the mechanical mode and with the anti-Stokes mode, the anti-Stokes 
mode is not entangled with the mechanical mode, but the whole system is 
in a fully tripartite-entangled state for a wide parameter regime. What is 
important is that in the present cavity scheme, such a parameter regime is 
much easier to achieve with respect to that of the free-space case. 
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Two driven modes and a 
movable mirror 

8.1 Introduction 

A bichromatic driving of a cavity has been experimentally realized by |57j . 
but their principal aim was to create an optical spring effect useful to trap 
and cool a heavy (m — lg) movable mirror. They theoretically predicted also 
the presence of optical entanglement between the output fields and intracavity 
light-mirror entanglement [12] • However, they computed the output optical 
entanglement only at the central lasers frequencies. 

Here we will use the theory of filters functions, developed in the previous 
chapter, to select physical output modes centered at different frequencies and 
with different bandwidths. Moreover the theory of filters will give us the 
possibility to calculate also the optomechanical entanglement of the output 
fields with the mirror, in a fully consistent way. 

In the previous chapter, we have seen that simply driving a cavity with 
a single laser we can get good optomechanical entanglement and also optical 
entanglement between the Stokes and anti-Stokes sidebands. That setup has 
the advantage to be very simple from an experimental point of view. However 
in this chapter we will see that, if we drive a cavity with two different lasers, a 
cooling laser (A a = u> m ) and a heating laser (As = -w m ), we can get better 
results. 

It turns out that optomechanical entanglement is stronger for the output 
mode associated with the heating laser. This is predicted in the rotating 
wave approximation, where entanglement is achieved for a field rotating with 
a frequency opposite to that of mirror. This means that, for every cavity 
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detuning, only the Stokes sideband is entangled with the mirror. Therefore, 
we can imagine that if the Stokes sidebend is on resonance (A# = —u> m ), then 
we can get a better optomechanical entanglement. This is true indeed, but 
we have already seen that at negative detunings the system is very unstable, 
since negative detuned lasers tends to drive the mirror. The idea is then to 
use another driving laser, with equal power and opposite detuning = uj m , 
to balance the heating/cooling rate and make the system stable. 

In this regime, both Stokes sidebands of the two lasers are entangled with 
the mirror and the entanglement of the mode driven by the heating laser is 
stronger, because of the cavity resonance. Moreover we get good results also 
for the optical entanglement of the output fields, where choosing the detection 
frequencies resonant with the cavity eigen-modes and using narrow detection 
bandwidths, we obtain a significant entanglement also at room temperature. 

8.2 Quantum Langevin equations 

We consider a Fabry-Perot cavity in which one of the end mirror is a microme- 
chanical oscillator. Two resonant modes of the cavity, oj c a and uj c b, are driven 
by two lasers with detuned frequencies ojqa and u>qb and powers Pa and Pb 
respectively. The Hamiltonian of the system can be written as, 

H = hw cA a^a + hw cB h% + ^huj m (p 2 + q 2 ) (8.1) 

-h(G A a* a + Gob tfb)q 

+ih[E A {a)e- iw » At - ae iuJ0At ) + E B {tf e^™ 1 - be*" 03 *)]. 

In the first line there is the energy of the optical modes and the mirror, which 
is modeled like a harmonic oscillator of frequency w m . The following commu- 
tation rules are satisfied, [a,aJ] = [b,tf] = 1 and [q,p] = i. The second line 
describes the interaction between the mirror and the light due to radiation 
pressure, with coupling constants Gq x = ^Jh/muj m uj cx /L, where m is the ef- 
fective mass of the mechanical oscillator, L is the length of the cavity, and 
x = A, B. The last line gives the contribution of the two driving lasers, in 
which \E X \ = y / 2P x n/huJo x , where k is the decay rate of the cavity, supposed 
to be the same for the two modes. In this treatment, scattering of photons 
of the driven modes into other cavity modes is neglected; this is acceptable if 
the mechanical frequency is much smaller then the free spectral range of the 
cavity, that is u; m « c/L. 

In interaction picture with respect to hwQA^a + fkvoB^b, we can derive a 
set of nonlinear QLE in which optical and mechanical noises are taken into 
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account, 

q = uj m p, (8.2a) 

p = -uj m q-'y m p + GoAa i a + GoBb' t b + £,, (8.2b) 

a = -[n + i(A A-Go A q)]a + E A + V2^a in , (8.2c) 

b = -[K + i(A B-Go B q)}b + E B + V2^b in . (8.2d) 

where Aq x = uo cx — ujq x are the detunings of the two lasers. The mechanical 
mode is affected by a viscous force with damping rate 7 m and by a Brownian 
stochastic force £(t), with zero mean value and the same correlations as given 
in £LU). 

The optical modes amplitudes instead decay at the rate n and are disturbed 
by the vacuum radiation input noises a m (t) and b in (t), characterized by zero 
mean values and the following correlations 

{a in {t)a in \t')) = {b in (t)b in \t')) = S(t - t'). (8.3) 



Setting the time derivatives to zero and solving for the mean values a 
(a), b s = (b), q s = (q), p s = (p), we get 

Ea 



E B 



(8.4) 
(8.5) 



G A\a s \ 2 + G 0B \b s \ 2 , . 

q s = , (8.6) 

Ps = 0. (8.7) 

where effective detunings A x = Aq x — (GQ A \a s \ 2 + GQ B \b s \ 2 ) / uj m have been 
defined, so that (j8.4H8.5p is actually a nonlinear system, whose solutions gives 
the stationary amplitudes a s and b s . 

For every operator O, we can define the deviation with respect to the 
steady state as 50 = O — (O). If we apply this transformation to equations 
(|8.2p and consider only first order terms in fluctuations, we obtain a set of 
linearized QLE for the deviation operators, 



5q = to m 5p, (8.8) 
5p = -uj m 5q - ~i m bp + G A \a 3 5<J + a*5aj 

+G 0B (b,6tf + b* s 5b) + e, 

5a = — (k + iA A )Sa + iGqa o-s dq + V^ko 111 , 
5b = -(K + iA B )5b + iG Ab s 5q + V2Kb in . 
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To simplify the following steps, it is convenient to define the effective coupling 
constants Ga = GqaclsV^, Gb = GosbsV^ and to choose a reference frame 
for the cavity modes so that a s = \a s \, and b s = |6 s |o 

It is better to work with field quadratures, which are defined for the mode 
"A as 5X A = (5a + 5a))/\f2, 5Y A = (Sa - 5a))/iy/2 and similarly for all 
the other bosonic operators. Now we can finally rewrite equations (18. 8p in 
matrix form as ii(t) = Au(t) + n(t), where u = (5q, dp, SXa, SYa, SXb, SYb) t , 
n = (0, £, 5Xf, 5YX 1 , 5Xf, 5Y l B n ) T , and 
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— K 


Ab 




\ 


Gb 











-A B 


— K 


J 



The differential equation can be formally integrated and gives the solution 



u(t) = M(t)u(0) + f dsM(t- s)n(s), (8.10) 
where M(t) = exp(At). 



8.3 Stability of the steady state 

The steady state exists and is stable if all the eigenvalues of the A matrix have 
negative real parts, so that M(oo) = 0. The characteristic polynomial of A is 



1 The linearized dynamics of this system is completely equivalent to that of the cavity with 
a membrane inside. In fact, the only qualitative difference between the two Hamiltonians 
(|6.69|) and (|8.2|) is in the sign of the radiation pressure term of the mode B. So, for the 
"membrane configuration" , we could repeat the same procedure obtaining the same equations 
but with the substitution Gob — > —Gob. However, by choosing the reference frame of the 
mode B such that b s = — |6 S |, the linearized QLE would be exactly equal to (|8.8p . 
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P(A) = A 6 + ciA 5 + c 2 A 4 + c 3 A 3 + c 4 A 2 + c 5 A + c 6 , where 

cl = 7m + 4K, (8.11) 

c2 = A 2 A + A 2 B + A-/ m K + 6k 2 + u4, 

c3 = 7m(A 2 4 + A| + 6/ t 2 ) + 2/ t [A 2 4 + A| + 2(k 2 + ^)], 

c4 = k 4 + 2 7m K(A| + 2k 2 ) +6k 2 w^ + A|(k 2 + u4) + 

Ai(A| + 2 lm n + k 2 + u 2 rn ) - to m {G 2 A A A + G 2 B A B ), 
c5 = 7 m (A 2 4 + K 2 )(A 2 3 + K 2 ) + 2^(Ai + A 2 + 2K 2 ) 

-2KUJ m (G 2 A A A + G 2 B A B ), 
c6 = uj 2 n (A 2 A + K 2 )(A 2 B + K 2 )-u; m [G 2 B A B (A 2 A + K 2 ) 

+G\A A {A 2 B + k 2 )]. 

We may ask if there exists a particular relation between the parameters, 
such that all the terms containing Ga and Gb in the previous coefficients 
(18. lip cancel out. In this case, the eigenvalues of A would be independent 
from the power of the two lasers and, in particular, these would be the same 
eigenvalues of the system with switched off lasers. The stability would be 
necessarily guaranteed. 

Directly from equations (|8.1ip . comes out that the eigenvalues of A are 
independent form Ga and Gb if and only if 

\G a \ = \G B \ = G, (8.12a) 
A A = -A B = A. (8.12b) 

If the these constraints are satisfied, the system is stable for any values 
of G and A. If instead (I8.12p are not satisfied, we have to apply the Routh- 
Hurwitz criterion. The inequalities that come out from this criterion are too 
involved to be reproduced here. 

The particular parameters configuration (|8.12p has a simple physical in- 
terpretation. Without loss of generality we can suppose A > 0, so that the 
system ()8.12p represents a perfect balance between a cooling laser (A^ > 0) 
and a heating laser A# < 0. 

This is actually an equilibrium regime very near the stability threshold, in 
fact it is dangerously broken as soon as the heating rate becomes higher than 
the cooling one. However we can easily avoid this problem by choosing the 
power of the heating laser slightly smaller than the cooling one. In this way 
one remains with certainty in the stable regime, also in the presence of power 
fluctuations of the driving lasers. 
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8.4 Entanglement measurement simulations using 
filter functions 

The two principal parameters which characterize a measurement of an out- 
put mode are the detection frequency Q and the bandwidth 7 which is also 
connected with the measuring time r ~ 1 /j. Such a measurement can be mod- 
elled using the theory given in the previous chapter, where the output field is 
filtered with a normalized function g(t) oscillating with a mean frequency 0, 
and with a bandwidth 7. 

In our particular system, we want to investigate two output modes orig- 
inating form two different cavity modes, therefore we do not need to choose 
orthogonal filter functions like that used for the single driving laser. 

The simplest choice is then to filter the two output modes a™ 1 (t) , a™ 1 (t) 
with two normalized damped plane waves with frequencies SIa,b and the same 
decay rate 7, 



9A(t) 
9B(t) 



2 ie - QAi ~' lt 0{t) , 



(8.13) 
(8.14) 



where 6{t) is the Heavyside step function. 
The two filtered modes are then 




(8.15) 



(8.16) 



which Fourier transformed gives 




(8.17) 



(8.18) 



where gA,B(w) are the Fourier transforms of the filter functions 




(8.19) 



7 + i(£Ia - u) ' 



(8.20) 



7 + i(Vt B - u) ' 



Now, in analogy with the previous chapter, we define: 
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a quadrature vector containing the mirror operators and the quadratures 
of the filtered modes, 

iT 



the corresponding correlation matrix, 

1 



V 



out I 



out 



(*)) 



.21) 



3.22) 



W = jW(t) + tT(t)' 

the transformation matrix containing the filter functions, 



T(t) = 



(8.23) 

• the matrix associated to the noise operators, which differently form the 
single mode case is diagonal because the input noise of the two optical 
modes are uncorrelated, 

/ o oooo\ 



( S(t) 
























S(t) 
























V2KRjeg A (t) 


-\/2K,lmg A {t) 


















^/2KlmgA (t) 


V2^Reg A (t) 
























\/2nRegB 


(t) 


— V2KlmgB(t) 




V o 











V / 2~Klm g b 


(t) 


V2KReg B (t) 


) 



D(cj) 








V o 



7m W 
(dm 





coth 







2k B T 





K 

k 

k 

k/ 



3.24) 



Repeating the same procedure used in the previous chapter, we find that V out 
is stationary in time, ad it is given by the same simple formula 



V 



out 



dwT{ui) 



M{u) + 



P< 



out 



2k 



D(u) 



P< 



out 



2k 



125) 



where P ou t = Diag[0, 0, 1, 1, 1, 1] is the projector onto the optical quadratures, 
T(lo) is the Fourier transform of (18.230 . and M(oj) if the Fourier transform of 
M{t) defined in rfHTUD . 

M(u) = (iuo + Ay 1 . (8.26) 

From the correlation matrix V out we can extract all the informations about 
the steady state of the system. In particular, choosing realistic parameters, 
we can compute the entanglement between the various lxl bipartite state 
obtained tracing out one of the optical output fields or the mechanical mode 
of the mirror. 
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8.4.1 Cooling mode-mirror optomechanical entanglement 

First of all we consider the entanglement between the cooling mode (A) and 
the mirror. The situation is very similar to the case that we have already 
considered with only one driving cooling laser. We choose the parameters very 
similar to those of the previous chapter, i.e., an oscillator with uj m /2TT = 10 
MHz, Q = 10 5 , mass m = 50 ng, a cavity of length L = 1 mm with finesse 
F = 2 x 10 4 , driven by two lasers with input powers Pa = 15 mW, Pb = 13 
mW, detunings = u m , A B = —uj m and wavelengths both around 810 nm. 
We have again assumed a reservoir temperature for the mirror T = 0.4 K, 
corresponding to n ~ 833. 

We are in the stable regime defined in (18.120 . but the power of the heating 
laser Pb is chosen slightly smaller than Pa to be far from instabilities. In fact 
all the eigenvalues Aj of A have negative real parts and the maximum of them 
is well far away from zero: max{Re(Aj)} ~ —2 x 10 5 Hz. 

In (Fig l8.ip a plot of the entanglement logarithmic negativity is given as a 
function of the measured frequency, for three different detection bandwidths 
determined by the parameter e = io m h- 

E;v 

0.20 r 




Figura 8.1: Logarithmic negativity Etf of the CV bipartite system formed 
by the mechanical mode and the cooling output beam versus the measured 
frequency Qa/^tu a t three different values of the inverse detection bandwidth 
e = w m /7. As expected only the Stokes sideband CI a = — w m is entangled 
with the mirror but, since it is off resonance, we obtain a small entanglement. 
The other parameters are uj m /27r = 10 MHz, Q = 10 5 , m = 50 ng, L = 1 mm, 
F = 2 x 10 4 , P A = 15 mW, P B = 13 mW, A A = uj m , A B = -u m , A = 810 
nm, T = 0.4 K. 
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The behavior of the optomechanical entanglement with respect to temper- 
ature is given in (Fig J8.2j) . 




T( K) 



Figura 8.2: Logarithmic negativity Ej^ of the CV bipartite system formed 
by the mechanical mode and the Stokes sideband (Qa = -^m) of the cooling 
laser versus the reservoir temperature T, at three different values of the inverse 
detection bandwidth e = ui m /j. The other system parameters are the same 
as in Fig. 18.11 



If we compare (Fig l8.ip with the monochromatic case (Fig l7.5p . we observe 
that the addition of a heating driving laser does not improve the situation, but 
it introduces more noise, decreasing the entanglement of the cooling mode. 



8.4.2 Heating mode-mirror optomechanical entanglement 



The advantage of the bichromatic driving is instead significant if we consider 
the entanglement of the heating laser beam (B) with the mirror. The entangled 
sideband, as predicted in the RWA, is again the Stokes one, but this time it is 
resonant with the cavity, in fact A = — uj m implies u> c a = uoa — w m = ui stokes- 

In (Fig l8.ip a plot of the entanglement logarithmic negativity is given as a 
function of the measured frequency, for three different detection bandwidths. 
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Figura 8.3: Logarithmic negativity Etf of the CV bipartite system formed 
by the mechanical mode and the heating output beam versus the measured 
frequency f2#/u; m at three different values of the inverse detection bandwidth 
e = w m /7- The Stokes sideband Qb = -w ra , which is on resonance with the 
cavity, is very entangled with the mirror. The other system parameters are 
the same in Fig. 18.11 

The behavior of the optomechanical entanglement with respect to temper- 
ature is given in (Fig J8.4j) . By choosing appropriate detection bandwidths, 
a significant amount of entanglement persists even above liquid He tempera- 
tures. 




Figura 8.4: Logarithmic negativity Ej^ of the CV bipartite system formed 
by the mechanical mode and the Stokes sideband (f2# = — u> m ) of the heating 
laser versus the reservoir temperature T, at three different values of the inverse 
detection bandwidth e = uj m /l- The other system parameters are the same 
as in Fig. 18.11 
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8.4.3 Heating mode-cooling mode optical entanglement 

Let us consider now the purely optical entanglement between the two output 
light beams. We optimize the system parameters, by choosing a cavity with 
a higher finesse J- = 8 x 10 5 and a lighter movable mirror m = 10 ng. We 
increase also the power of the driving lasers to Pa = 75 mW and Pg = 65 mW. 
We are also in this case near the stable situation of (I8.12p . with Pb (heating 
laser) slightly smaller than Pa (cooling laser). 

Stability is confirmed by the eigenvalues of A, all having negative real parts 
with max{Re(Aj)} ~ -5 x 10 5 . 

The RWA suggests us to consider, also in this case, counter-rotating modes 
to achieve a better entanglement. In fact it comes out that the best choice is to 
consider the anti-Stokes sideband of the cooling beam and the Stokes sideband 
of the heating beam, that is Qa = and SIb = —^m- We observe that also 
in this case the entangled sidebands are both resonant with the cavity. 

In (Fig l8.5p . the logarithmic entanglement negativity is plotted versus the 
measured central frequency Qa/^tti of the cooling mode, while the detected 
frequency of the heating mode is centered on the Stokes sideband Qb = — w m . 




Figura 8.5: Logarithmic negativity Ej^ of the CV bipartite system formed 
by the cooling and the heating output beams versus the measured frequency 
Q,A/^m of the cooling beam at three different values of the inverse detection 
bandwidth e = uj m /l- The detected frequency of the heating beam is centered 
on the Stokes sideband. The system parameters are ui m /27r = 10 MHz, Q = 
10 5 , m = 10 ng, L = 1 mm, F = 8 x 10 5 , P A = 75 mW, P B = 65 mW, 
A A = uj m , A B = -u m , X = 810 nm, T = 0.4 K. 

Differently form the optomechanical entanglement, the logarithmic nega- 
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tivity associated to the two optical output beams, always increases for decreas- 
ing detection bandwidths. However the more the bandwidths are narrow, the 
more the measured frequencies must be centered with high precision and this 
imposes a realistic upper limit for the parameter e. 

By measuring narrow bandwidths, thermal noise can be reduced, giving 
a significant amount of optical entanglement even at room temperature (see 
FigEH). 




Figura 8.6: Logarithmic negativity Ej^ of the CV bipartite system formed 
by the anti-Stokes sideband of the cooling mode (Qa = ^m) an d the Stokes 
sideband = — u m ) of the heating laser versus the reservoir temperature 
T, at three different values of the inverse detection bandwidth e = uj m /^. The 
other system parameters are the same as in Fig. 18.51 

8.4.4 Tripartite entanglement 

In analogy with the monochromatic configuration, we check the presence of 
tripartite entanglement between the Stokes sideband of the heating mode, the 
anti-Stokes sideband of the cooling mode and the mirror. 

With the same procedure used in the previous chapter, we consider the 
correlation matrices after partial transposition with respect to one of the three 
sub-systems and then we check if the Heisenberg condition (I2.22|) is violated. 
In Fig. 18.71 the minimum eigenvalue of the three test matrices is plotted as a 
function of the inverse bandwidth e. Also in this bichromatic setup, we find 
a large region in which the system is in a fully tripartite-entangled state (all 
the three test matrices are not positive definite). 
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Figura 8.7: (Color online) Analysis of tripartite entanglement. The minimum 
eigenvalues after partial transposition with respect to the Stokes sideband 
(Ob = — Wm) of the heating mode (blue line), anti-Stokes sideband (£Ia = u m ) 
of the cooling mode (green line) and the mechanical mode (red line) are plotted 
versus the inverse bandwidth e = u m /j. The other parameters are the same as 
in Fig. 18.11 These eigenvalues are all negative in the studied interval, showing 
that one has fully tripartite-entanglement. 
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Conclusions 



In the first part of the thesis some aspects of continuous variable quantum 
information has been investigated. 

After a general introduction on the phase space representation of Gaus- 
sian states, we have studied the quantum mechanical phenomenon of entan- 
glement restricted to CV systems. The most relevant separability criteria and 
entanglement measures have been given. 

In chapter^ the Braunstein-Kimble protocol for the teleportation of Gaus- 
sian states has been described in the Heisenberg picture, while in chapter [5j we 
have presented some original results on the optimization of the teleportation 
fidelity by local Gaussian operations. We have shown that, for a given shared 
entanglement, the maximum fidelity of teleportation is bounded below and 
above by simple expressions depending upon the lowest PT symplectic eigen- 
value v only (see Eq. (|5.ip ). We have seen that these bounds are quite tight 
and that the upper bound of the fidelity is reached if and only if Alice and Bob 
share a symmetric entangled state. We have also determined the general form 
of the CM of Alice and Bob state after the optimization procedure. Then we 
have restricted to local TGCP maps and shown that the optimal TGCP map 
is composed by a local symplectic map, eventually followed by an attenuation 
either on Alice or on Bob mode. Finally we have shown how the corresponding 
value of the maximum fidelity J- op t can be derived from the knowledge of the 
symplectic invariants of the initial CV entangled state shared by Alice and 
Bob. The optimization by generic GCP maps (i.e., including also Gaussian 
measurements on ancillas) or by non-Gaussian local operations are still open 
questions. 

The second part of the thesis is all focused on optomechanical quantum 
systems. 

In chapter [6l we have derived the quantum Langevin equations related to 
optical cavities (input-output theory) and to mechanical systems (quantum 
Brownian motion). The radiation pressure potential between a mechanical 
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resonator and an optical cavity mode has been also derived through different 
approaches. 

In chapter [7] we have considered an optical cavity mode coupled with an 
oscillating mirror by radiation pressure. We have studied in details the light- 
mirror entanglement and its relationship with the laser-cooling of the mechan- 
ical mode. We have developed a consistent theory to define filtered output 
modes, with given central frequencies and detection bandwidths. We have 
used this theory to show that the Stokes and anti-Stokes sidebands (photons 
scattered by the mirror) are optically entangled, and that the Stokes one is 
also robustly entangled with the mirror vibrational mode. We have seen also 
that the two sidebands and the mirror mode form a fully tripartite-entangled 
state. 

In chapter [H] we have extended the previous model by using two driving 
lasers. We have observed that if the two lasers have equal powers and opposite 
detunings, then a balance between the cooling and the heating rate can be 
reached, giving the possibility to achieve interesting stable regimes. In this 
configuration, the Stokes sideband of the heating laser can be chosen resonant 
with the cavity and so it can reach a high amount of entanglement with the 
mirror. Moreover, two output modes associated to counter rotating sidebands 
(e.g. Stokes of the heating laser and anti-Stokes of the cooling laser) can 
achieve a significant optical entanglement even at room temperature, while at 
low temperature they form, together with the mirror vibrational mode, a fully 
tripartite-entangled state. 
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questo pezzo di storia ethmero? Solo l'Eterno, l'Eterno che ama l'emmero. In 
Lui spero. 
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